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DIGROUP  DATA  REDUCTION  TECHNIQUES 


1.  INTRODUCTION  AND  SUMMARY  OF  RESULTS 

The  purpose  of  this  study  Is  to  deterrlne  the  best  way  to 
conpress  the  output  of  tv;o  T1  PCM  carrier  terminals  down  Into  the 
bit  stream  normally  required  by  one  T1  terminal*.  The  equipment 
which  performs  this  2:1  compression  Is  herein  called  a Dlirroup 
Data  Reduction  (DDR)  system.  The  use  of  this  system  can  save 
money  because  It  doubles  the  number  of  speech  channels  that  can 
be  placed  on  a T1  digital  transmission  line  and  thereby  lowers 
the  per  channel  transmission  cost.  Even  when  fully  loaded  there 
Is  no  noticeable  degradation  to  the  user. 

A block  diagram  of  the  DDR  system  and  how  It  Is  used  Is 
shown  In  Figure  1.  Only  one  direction  of  transmission  Is  shown. 

A complete  DDR  terminal  consists  of  a DDR  encoder  and  a DDR 
decoder.  T1  carrier  systems  normally  require  transmission  over 
one  1.5^^  mB  line  In  each  direction.  By  uslne  the  DDR  terminal 
two  complete  T1  carrier  systems  require  a single  m.B  line 

In  each  direction.  As  shown  In  the  figure  the  DDR  system  can 
accommodate  up  to  six  dedicated  data  channels  on  which  no 
compression  is  possible  because  the  nature  of  the  data  Is  un- 
restricted. These  data  channels  are  obtained  by  designating  up 
to  six  of  the  T1  speech  channels  as  data  ports.  These  data  ports 
are  designed  to  accept  digital  data.  They  are  routed  straight 
through  the  DDR  system  without  compression.  Sv/ltches  on  the  DDR 
terminal  designate  which  of  the  T1  speech  channels  are  to  be 
regarded  as  data  channels  and  omitted  for  the  data  comoresslon 
al gorlth.m.s  designed  for  soeech. 

*A  T1  carrier  system  carries  2*i  two-way  speech  communication 
channels.  It  Is  widely  used  In  military  and  civilian  teleohor.e 
networks. 
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Chapter  2 of  this  report  suranarlzes  four  previously 
proposed  data  conrresslon  systenis  that  have  been  considered  In 
the  past  for  applications  like  this.  Of  these  four  only  one 
seemed  promising  - the  Speech  Predictive  Encoding  Ccmunlcat Ion 
(SPEC)  system  currently  being  tested  by  COMSAT.  In  addition 
to  SPEC  we  determined  that  two  other  promising  possibilities 
existed  - NIC/TASI,  a marriage  of  Nearly  Instantaneous  Companding 
(NIC)  with  Time  Assignment  Speech  Interpolation  (TASI)  and 
ADPCM/TASI,  a marriage  of  Adaptive  Differential  Pulse  Code 
Modulation  (ADPCM)  with  TASI.  These  three  system.s  - SPEC,  NIC/ 
TASI  and  ADPCM/TASI  - were  considered  to  be  candidates  for  the 
DIP.  application  and  are  given  further  study  In  Chapter  IV.  The 
last  section  of  that  chapter  contains  our  reasons  for  choosing 
ADPCM/TASI  as  the  system  best  suited  for  the  DDP.  application. 

The  COMSAT  SPEC  system  Is  reasonably  well  documented  [1,2]. 
It  Is  currently  undergoing  tests  on  satellite  links.  The  concept 
of  ADPCM  Is  well  known  [3f**]  tut  Its  use  with  TASI  has  not  teen 
previously  explored  In  any  detail.  The  TASI  system  Itjelf  has 
long  teen  used  on  undersea  analog  cable  transmission  systems  [5,6 
The  NIC  system  Is  the  newest  speech  data  compression  teohnlcue 
considered  herein.  It  was  presented  at  a conference  In  1?75  [’’[ 
as  a possibility  fcr  mobile  radio  communications  but  the  appli- 
cation to  cable  car;  ler  systems  was  not  evident  until  Septe.mber 
1976  [8]  - only  about  two  months  before  the  termlr.atlcn  of  this 
study.  All  of  the  t.hree  candidate  compression  systems  take 
advantage  of  two  characteristics  of  speech:  (1)  the  user  of  a 
tvra-way  comnur.lcat lor.  system  talks  only  about  33'  cf  time  the 
system  Is  off-hook  (^.e. , the  duty  factor  Is  .39)  and  (2)  speech 
sample  values  are  slowly  varying  comoared  to  the  PCM  sa.mpll.ng  rat 
a.nd  are  therefore  somewhat  predictable. 
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since  all  three  candidate  systen*  take  advantage  of  the 
average  duty  factor  of  speech,  which  Is  about  .38.  they  must 


adapt  to  those  times  when  the  Instantaneous  duty  factor  Is 

greater  than  .38;  there  are  times  when  the  system  Is 

heavily  loaded  and  more  than  38!  of  the  users  In  one  direction 

of  transmission  are  talking.  The  SPEC  system  adapts  to  this  by 

transmitting  only  the  most  unpredictable  speech  sam.ples  during 

high  activity.  Those  samples  which  are  frozen  out  are  predicted 

by  the  decoder  to  be  Identical  to  the  most  recently  transmitted 

value.  The  ADPCM/TASI  and  NIC/TASI  systems  adapt  to  high  activity 

by  reducing  the  number  of  bits  used  to  encode  each  sample  value. 

This  technique  Is  used  In  Europe  on  some  PCM  links. 

The  principle  tool  used  In  this  study  Is  computer  simulation. 

^ Analyses  of  SPEC,  ADPCM,  TASl  and,  to  a lesser  degree  NIC,  are 

I 

available  In  the  publlsned  literature  [1-8].  What  we  needed  to 
determine  was  how  these  systems  sound  in  a system  designed  for 
the  DDR  application.  Computer  simulation  allowed  us  to  program- 
matically vary  a large  number  of  the  syste.m  parameters  - such  as 
system  loading,  duty  factor.  Input  signal  level,  number  of 
dedicated  data  channels,  etc.  The  simulations  were  done  or.  an 
AOT/30  Interactive  computer.  Telephone  quality  speech  was 
sampled,  ceded  and  stored  on  a digital  disk.  Signal  processing 
software  which  simulated  various  DDR  systems  was  written  (In 
Fortran  IV)  and  the  speech  was  passed  through  these  programs. 

The  processed  speech  could  then  be  listened  to  and  recorded  on 
audio  tape.  Nc  formal  subjective  tests  were  made  but  an  audio 
demonstration  tape  was  r.ade  and  forms  a part  of  this  report. 

This  audio  tape  Is  not  furnished  with  every  report  but  a copy 


of  the  tape  can  be  obtained  If  requested  before  November  31 » 1977 
- see  Chapter  6 for  ordering  Information  and  a description 
of  the  tape.  The  computer  programs  also  furnished  a number  of 
objective  performance  measures  such  as  slgnal-to-noise  ratios. 

The  computer  simulations  were  used  to  compare  three  CLR  systems: 

(1)  A lightly  loaded  system,  with  only  16  off-hook  users.  The 
system  performs  test  under  this  condition  because  the  speech 
compression  algorlthm.s  are  not  stressed.  This  system  Is 
called  the  l6/0  system. 

(2)  A fully  loaded  system,  with  k6  off-hook  telephone  users  and 
two  dedicated  data  channels.  This  system  Is  called  the  ^6/2 
system. 

(3)  A fully  loaded  system  with  1)2  off -hook  telephone  users  and 
six  dedicated  data  channels.  This  is  called  the  1(2/6 
system.  The  speech  compression  algorithms  are  highly 
stressed  In  this  configuration.  It  represents  a worst  case 
condition  - the  speech  compression  Is  2.3  to  1. 

On  the  basis  of  these  computer  simulations  and  other  con- 
siderations we  feel  that  the  ADPC.M/TASI  system  Is  best  for  the 
ODR  application.  This  conclusion  reflects  the  fact  that  the 
ADPCM  Is  tetter  suited  for  compressing  speech  signals  because  It 
Is  more  "tailored"  to  the  statistics  of  speech  than  the  other 
two  candidate  systems.  The  ADPCM  system  Is  not  likely  tc 
operate  well  with  non-speech  signals  such  as  voice  band  data. 

Ve  have  assumed  that  all  data  Is  presented  tc  the  data  ports  In 
digital  form  and  that  voice  band  data  Is  not  present  on  the 
speech  channels. 
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2.  REVIEW  OP  PREVIOUS  DIGITAL  SPEECH  COMPRESSION  SYSTEMS 

Currently  there  Is  a lot  of  activity  In  the  teleconununlca- 
tlons  Industry  In  the  design  of  systems  for  digital  speech 
compression.  Very  complicated  systems  for  this  application  can 
be  constructed  for  reasonable  cost  using  LSI. 

Table  1 Is  a compilation  of  four  systems  that  have  been 
previously  proposed  for  multichannel  bit  rate  compression  of 
I digitally  encoded  speech.  The  last  system  In  the  table,  SPEC, 

Is  being  actively  considered  for  applications  In  the  telecommuni- 
cations Industry  now  and  Is  a candidate  for  the  DDR  application. 
Following  this  table  are  brief  outline  descriptions  of  the  systems 
listed  In  the  table. 
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2.1  Vocoplexer  or  ODVP  (Optimum  Digital  Voice  Processor) 
References:  [12,13] 


NOTES 

1*  Research  sponsored  by  ECOM  and  RADC. 

2.  Digital  data  over  troposcatter . 

3.  Constructed  by  Bell  Aerospace. 

ll.  Uses  pauses  In  speech  — 25*  duty  factor. 

5.  Codes  output  of  6 bit  PCM  to  give  60S  reduction  — *i8  speech 
channel  PCM. 

6.  Run  length  coding  starting  with  most  slg  bit  of  each  channel. 
Probably  very  sensitive  to  errors  — cne  error  can  cause 
errors  In  the  rest  of  the  frame  (one  bit  error  can  cause 
catastrophic  failure  In  all  channels).  Transmits  sign  bit 

as  Is. 

7.  Experiment  was  performed  (field  test). 

8.  Effect  of  buffer  overflow  Is  to  reduce  system  to  5 bit  PCM 
operation  — graceful  degradation. 

9.  Hardware  model  tested  with  18  channel  fully  Implemented  plus 
simulator  to  simulate  .25  activity  factor.  Also  one  2.ii  kB 
data  channel . 

10.  Designed  to  operate  without  compander. 

11.  Effect  of  errors  not  considered  In  test  model. 


2.2  PCM/TASI 

References:  [10,11] 


NOTES 

1.  08  channel  6 bit  PCM  over  TROPO  (AN/MRC-98). 

2.  Compresses  08  channels  dovm  to  19  ♦ 1 channel  overhead. 

3.  Average  length  of  freezeout  Is  32  m sec. 

0,  Each  125  usee  frame  has  19  6 bit  channel  blocks  » 8 bits 
for  channel  signaling  which  specifies  channel  assignment. 

1 multi  frame  » 00  frames. 

.•.  There  are  00  x 8 ■ 320  bits  for  channel  assignment. 

« l6  channel  assignment  bits  per  channel  each  multlframe. 
Use  ^5^ehan  bits  redundancy. 

5.  BCH  code  with  7 Info  + 8 redundancy  + 1 parity  check  bit  for 
7 Info  bits.  Hanning  distance  of  6 .•.  3 bit  errors  can  be 
corrected.  When  error  Is  detected  channel  Is  squelched  — 
this  means  that  channel  assignment  Is  not  updated  — previous 
channel  asslgr.ment  Is  used.  This  Is  like  SPEC. 

6.  Channel  assignment  updated  every  5 msec. 


2.3  Direct  Entropy  Coding  of  PCM  Levels 

Reference  : [9] 

NOTES 

1.  Assign  a slnpliried  entropy  code  to  PCM  channels 

transmit  0 for  Inactive  channels  1 called  a 

transmit  1+8  bit  code  for  active  channels  i 1-9  code 

2.  Achieve  a rate  of  l4.ll  bits  per  8 bit  sample. 

3.  Other  entropy  codes  are  possible  — the  1-9  code  Is  simple. 
1).  l6  bit  framing  and  signaling  word. 

5.  Like  SPEC,  they  force  a zero  code  In  Inactive  channels. 

6.  Do  not  eliminate  predictive  redundancy  from  signal  — use 
only  activity  factor. 

7.  No  channel  signaling  information  Is  considered. 
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2.H  SPEC 


References : [1,2] 


NOTES 

1.  Converts  6*1  speech  channels  down  to  32  — 2.0U8  m3/sec. 


2.  Unlike  other  syster.s  It  takes  sone  advantage  of  predicta- 
bility of  speech. 

3.  Digital  vclce  switch  used  to  ell.-nlnate  noise  during  silent 
Intervals.  Only  active  speech  goes  through. 

t.  Uses  zero  order  predictor. 

5-  Assignment  recirculation  unit  rotates  priority  of  service 
to  Input  trunks. 

6.  6A  bit  sample  assignment  word  (SAW)  contains  a 1 for  each 
active  channel. 


7.  Sample  Is  transmitted  only  If  It  exceeds  previous  sample 

by  some  6 — 6 ■ 3 is  Im.perceptable  from  6-0.  6*3 

means  3 quantizing  levels. 

8.  Service  priority  circulated  over  all  trunks  to  give  unifcrr 
degradation  to  all. 


9.  Single  parity  check  on  SAW. 


10. 


11, 


Uses  A law  companding,  A • 87.6. 

Activity  factor  with  « - 0 is  33*1  q,  chanire 
Activity  factor  with  6 - 3 Is  2UXJ  cnange. 

S Is  called  predictor  aperture  or  threshold.  5 changes 
from  0 (low  loading)  to  11  (high  loading). 


12.  Busying  out  channels  la  not  used  — degradation  Is  spread 
uniformly  Into  all  speech  channels  — there  Is  no  priority 
system.  Data  channels  are  not  handled  differently  from  speech. 


13.  Spec  Is  equivalent  to  standard  8-blt  PCM  95*  of  the  time  — 

field  test  (must  vary  with  # and  type  of  VB  data  sets  on  line). 

1^.  Intelligible  crosstalk  cannot  occur. 
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3.  THE  TASI  PRINCIPLE 


TASI  (Time  Assignment  Speech  Interpolation)  Is  a technique 
which  allows  c transmission  channels  to  be  shared  among  n speakers 
(or  usjrs).  It  takes  advantage  of  the  fact  that  a user  of  a 
speech  communications  channel  talks,  on  the  average,  only  about 
381  of  the  time.  For  the  remaining  62J  of  the  time  the  user  Is 
either  listening  or  pausing.  Thus  the  user  needs  a transmission 
channel  or  path  only  38I  of  the  tlm.e.  TASI  takes  advantage  of 
this  38<  duty  factor  by  detecting  when  speech  Is  present  and 
connecting  the  speaker  to  a transmission  channel  only  when  he  Is 
speaking.  Extensive  studies  have  been  made  on  talker  statistics 
and  the  advantages  to  be  gained  using  TASI  [5,6]. 

The  degradation  usually  associated  with  TASI  systems  Is 
clipping  on  the  Initial  sounds  In  a speech  burst.  The  subjective 
effect  of  speech  clipping  Is  negligible  If  the  clipping  Is  less 
than  about  50  m sec  [l*t].  Thus  the  Job  of  the  TASI  system  Is  to 
determine  when  speech  Is  present  and  switch  the  speaker  to  a 
transmission  channel  within  50  m sec.  Some  fraction  of  the  bit  rate 
In  a digital  TASI  system  must  be  devoted  to  transmitting  this 
voice  switch  Information  from  encoder  to  decoder  so  that  the 
decoder  knows  which  speaker  Is  connected  to  each  of  the  trans- 
mission channels.  This  voice  switch  Information  is  contained  In 
what  we  call  the  "channel  assignment  message". 

In  the  DDR  application  one  frame  Is  193  bits  long  and 
lasts  for  125  usee.  At  least  once  every  50  m sec  (1400  frames) 
the  channel  assignment  message  must  Inform  the  decoder  of  which  of 
the  ^6  speech  channels  Is  on  and  should  be  assigned  to  a transmission 
channel.  A simple  *16  bit  word  can  do  this  by  transmitting  a 1 In 
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each  of  the  positions  If  that  voice  switch  Is  on.  Otherwise 
a zero  Is  transmitted.  This  requires  only  an  average  of  • 

.11  bits  per  frame.  In  the  systems  described  herein  we  devote 
two  bits  per  frame  for  this  voice  switch  Information  and  two 
bits  for  error  protection.  Using  a rate  j error  correcting  code 
for  error  protection,  which  Is  needed  because  errors  In  the 
assignment  word  are  very  disruptive,  ^ bits  per  frame  are  used 
for  the  channel  assignment  message.  This  allows  the  voice  switch 
Information  to  be  transmitted  once  every  23  frames  (2.9  msec). 

In  what  follows  this  bits  per  frame  Is  allocated  for  the  channel 
assignment  message  for  both  the  and  *12  speech  channel  versions 
of  the  DDR  system. 

3.1  Signaling  over  a Digital  TASI  System 

The  PCM  systems  which  are  the  Inputs  to  the  DDR  compressor 
send  one  signaling  bit  every  sixth  frame  or  1/6  bit  per  frame. 

For  *J6  Input  channels  this  Is  about  1/6  x **6=8  bits  per  frame. 

In  the  TASI  systems  simulated  and  studied  herein  we  have  assigned 
8 bits  of  each  193-blt  frame  to  be  signaling  bits.  The  signaling 
Information  Is  stripped  off  the  *l6  PCM  channels  and  transferred 
to  this  8-blt  TDM  signaling  "channel". 

In  this  method  the  DDR  system  Is  transparent  to  the  signaling 
Information.  In  a more  sophisticated  DDR  system  we  could  take 
advantage  of  the  fact  that  the  PCM  signaling  Information  has  a 
very  low  Information  rate.  The  maximum  signaling  rate  for  each 
speech  channel  Is  20  pulses  per  sec.  (actually  lOpps  Is  adequate 
for  almost  all  systems).  For  a 20pps  rate  each  dial  pulse  must 
be  switched  on  and  then  switched  off  — this  requires  **0  blts/sec 
for  *«6  channels:  **6  x **0=18*40  blts/sec  are  needed.  This  Is  only 
about  1 pulse  every  **  frames.  This  signaling  Information  must 
be  protected  against  transmission  errors.  Using  a 1/2  rate  error 
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correcting  cod*  would  require  1 bit  every  two  frames.  This 
could  be  accomplished  by  using  every  other  framing  bit  which 
In  PCM  systems  Is  reserved  for  Common  Channel  Interoffice 
Signaling  (CCIS).  The  T-carrler  systems  are  built  to  frame  on 
1 framing  bit  every  2 frames.  With  this  method  the  S-blt  per 
frame  now  required  for  signaling  could  be  used  for  speech 
transmission. 

In  an  even  more  sophisticated  arrangement  both  the  signaling 
and  voice  switch  Information  could  be  encoded  Into  one  6-blt 
word  per  frame.  This  word  could  then  be  transmitted  only  v;hen 
needed  - ^.e.,  only  when  voice  switch  or  signaling  Information 
Is  required.  In  whatever  signaling  scheme  Is  used  the  probability 
of  a signaling  error  should  be  much  lower  than  the  p(e)  provided 
on  typical  T1  lines.  The  telephone  company  uses  error  detection 
and  retransmission  for  all  of  Its  CCIS  channels.  Error  proba- 
bilities on  the  order  of  10"^°  are  achieved  on  signaling  channels 
of  the  tell  plant  of  the  Bell  System.  For  the  purpose  of  this 
study  we  have  assigned  8 bits  per  frame  to  a separate  slg.^allng 
channel  as  discussed  In  the  first  paragraph  of  this  section. 

3.2  The  Digital  TASI  Frame  ■» 

In  this  study  we  have  assumed  that  for  each  125  usee  frame 
the  voice  switch  Information  occupies  a i(-blt  channel  assignment 
channel;  the  signaling  occupies  8 bits;  and  the 
remaining  180  bits  are  devoted  to  data  and  speech  transmission. 

The  frame  tit  assignments  for  the  1*6/2  and  42/6  systems  are 
shown  In  Figure  2. 


3 . 3 Talker  Statistics 


The  purpose  of  this  section  Is  to  present  talker  statistics 
for  multiple  channel  carrier  systems  with  k6  and  1)2  off-hook 
speech  channels.  The  systems  described  herein  take  advantage 
of  the  fact  that  on  the  average  a one-way  channel,  when  In  use, 
contains  speech  only  38*  of  the  tine  [1].  The  other  62*  of 
the  time  the  user  Is  silent  - either  because  he  Is  listening  or 
pausing  during  a conversation. 

Assume  there  are  n Independent  talkers,  each  of  whom  has 
an  average  activity  factor  p.  At  any  particular  Instant  the 
probability  that  the  number  of  simultaneous  talkers  will  equal 
or  exceed  c Is  given  by  the  cumulative  binomial  distribution: 

"X  x-iuUyT  p' 

When  c and  n are  large  enough  B(c,n,p)  can  be  approximated 
by  the  cumulative  Gaussian  distribution  and  In  the  case  where  B 
Is  less  than  .50  and  n and  c are  large  B(c,n,p)  can  be  computed 
by  this  formula  [15]: 

B(c,n,p)  - Exp  [-  -^  (1  + y + .6y^)]  (1) 

where  y > 

0 • /np( 1-p) 
c • number  of  channel:, 
p • average  activity 
n ■ number  of  talkers 

This  formula  can  be  used  to  determine  what  percentage  of 
the  time  c channels  are  active  when  r.  channels  are  off-hook. 


t 
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For 
B(c ,n,p) 


systems  with  n”lt6  and 
Is  tabulated  In  Table 


1*2  off-hook  channels  the  function 
2 for  some  relevant  values  of  c. 


Table  2.  Tabulation  of  B(c,n,.38) 


c I n~l«6 I n*lt2 


ID* 

.67 

15* 

.66 

16* 

.65 

17 

.59 

.D3 

18 

.D9 

.32 

19 

.380 

.211 

20 

.272 

.130 

21 

.180 

.073 

22 

.110 

.038 

23 

.063 

.018 

2D 

.033 

.0079 

25 

.0159 

.0032 

26 

.0071 

.0012 

27 

.0029D 

.OOOD 

3D 

1 .8x10' 

•These  entries  obtained  from  B(c , n, . 38 ) -1-B( n-c , n , . 62 ) . 

For  systems  with  *<6  and  *<2  off-hook  channels  the  probability 
that  exactly  c channels  are  In  use  are  contained  In  Table  3 for 
various  values  of  c. 


Table  3*  Number  of  simultaneous  talkers  for 
system  with  1)6  off-hook  users. 


c Percentage  of  time  c users  are 

number  of  talking  simultaneously 


simultaneous 

talkers 

n«D6  off-hook  chan. 

[ n»D2  off-hook  chan 

16  or  less 

57 

17  or  less 

51 

68 

17 

10 

11 

18 

11 

11 

19 

11 

8.1 

20 

9.2 

5.7 

21 

7.0 

3.5 

22 

D.7 

2.0 

23 

3.0 

1.0 

2D 

1.7 

.D6 

25 

.88 

.20 

26  or  more 

.71 

.12 

This  table  shows,  for  example,  that.  In  a system  with  *16 


off-hook  speech  channels,  the  probability  Is  only  .03  that 
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23  people  will  be  talking  simultaneously.  The  systems  described 
herein  take  advantage  of  this. 

Bits  per  sample  on  TASI-type  systems 

Using  Table  3 we  car.  calculate  the  probability  that,  say, 

L bits  will  be  available  to  any  given  speech  sample.  Considera- 
tion of  voice  switch  requirements  has  led  to  the  conclusion  that 
for  the  U6/2  system  l6lt  bits  are  available  per  frame  for  the 
trans.mlsslon  of  the  speech  samples.  In  this  system  16  bits  are 
being  used  for  the  dedicated  data  lines  and  12  bits  are  used  to 
transmit  voice  switch  and  signaling  Information.  For  the  *12/6 
system  132  bits  are  available  for  speech  sa.mples,  are  used 
for  data  and  12  for  voice  switch  and  signaling  Information.  The 
probability  that  L bits  are  available  for  speech  samples  Is 
given  In  Table  ^4 . 


Table  4.  Probability  that  L bits  are  available 
for  each  speech  sample. 


L • bits  available 
for  speech  sample 

4)6/2  system 
(164)  bits  available) 

4)2/6  system 
(132  bits  available) 

8 

.82 

.57 

7 

.14)7 

.22 

.032 

.19 

5 

.018 

4)  or  less 

.0004) 

Examination  of  Table  i|  shows  that  for  the  i»2/6  system  6 or 
more  bits  are  available  for  each  speech  sample  98S  of  the  time. 

For  the  4)6/6  system  6 or  more  bits  are  available  99-9*  of  the  tine. 
3.4*  High  Activity  Factor  Statistics 

Although  the  systems  In  this  report  are  compared  assuming 
a 38t  activity  factor,  there  are  times  when  the  activity  factor 
exceeds  this  value  for  minutes  at  a time.  In  COMBAT 's  report 
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on  the  field  trial  of  their  SPEC  system  it  Is  shown  that  during 
one  three  minute  period  the  activity  factor  exceeded  U^%  (see 
Fig.  ll(b)  of  [2]).  There  were  several  other  shorter  periods  when 
the  activity  factor  was  about  *)5%.  In  mllllary  situations  It 
seems  likely  that  In  some  circumstances  the  activity  factor  In 
one  direction  or  another  might  exceed  U5!.  For  example.  If 
there  were  significant  military  activity  In  one  geographical 
region  then  telephone  circuits  coming  out  of  this  region  would 
tend  to  have  a higher  activity  factor  than  those  circuits 
directed  Into  the  region;  i-e..  If  this  geographical  region  Is 
generating  military  Information  then  users  In  the  region  may 
tend  to  be  more  active  talkers,  whereas  the  people  they  are 
conversing  with  would  tend  to  he  listeners. 

Although  this  Issue  Is  not  considered  further  In  this 
report  talker  statistics  have  been  computed  for  an  activity  factor 
of  45s  for  the  ii6  and  1)2  channel  systems  considered  herein. 

Table  5 shows  the  quantity  B(c,n,.l)5)  for  n»46  and  n=i)2.  Table  6 
shows  the  percentage  of  time  c users  are  talking  simultaneously 
for  these  two  systems. 


Table  5*  Tabulation  of  B(c,n,.il5) 


c 

n-i)6 

n*l)2 

19 

.535 

20 

.1)26 

21 

.513 

.312 

22 

.1)07 

.211 

23 

.300 

.1319 

21) 

.205 

.0758 

25 

.1298 

.0001 

26 

.0761) 

.01959 

27 

.01)17 

.00881 

28 

.02117 

.00365 

29 

.00996 

.001395 

30 

.OOU35 

.0001)02 

31 

.001765 

32 

.000661) 

19 


r 


c-nuBber  of 

simultaneous 

talkers 

Percentage  of  tl 
talking  slmultar 

n***6  off-hook  chan. 

me  c users  are 
leously 

n“**2  off-hook  chan. 

l8  or  less 

.1<65 

19 

.109 

20  or  less 

.1*87 

.688 

20 

.096 

.IIH 

21 

.106 

.101 

22 

.107 

.079 

23 

.095 

.056 

21* 

.0752 

.0357 

25 

.053'! 

.0205 

26 

.03‘*7 

.0108 

27 

.0205 

.00516 

28 

.0112 

.00225 

29 

.00561 

.0009 

30 

.00258 

31 

.00110 

Table  6.  User  statistics  for  the  high  activity  factor  of  .*<5. 
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4,  CANDIDATE  DDR  SYSTEMS 


In  our  Judgement  there  are  three  candidate  systems  which 
could  best  satisfy  the  task  of  combining  two  T1  24  channel 
terminals  Into  the  1.54mb  bit  stream  normally  occupied  by  one 
terminal.  These  are 

1.  SPEC  (Speech  Predictive  Encoding  Communication  system) 
proposed  by  CO.MSAT 

2.  NIC/TASI  (Nearly  Instantaneous  Companding,  a form  of 
adaptive  pulse  code  modulation  com.blned  with  Tine 
Assignment  Speech  Interpolation) 

3.  ADPCM/TASI  (Adaptive  Differential  Pulse  Code  Modulation 
combined  with  TASI) 

At  the  time  of  this  writing  COMSAT  has  published  a consid- 
erable amount  of  information  [1,2]  about  earlier  versions  of 
SPEC  although  recent  Impi  ovements  are  not  documented.  Nothing 
has  been  published  about  the  NIC/TASI  system  by  Bell  Telephone 
Laboratories  which  Is  apparently  performing  tests  on  such  a 
system  between  Hew  York  and  Boston  [16].  Two  articles  [7,8] 
discuss  NIC  tut  not  In  conjunction  with  TASI.  There  Is  consid- 
erable Information  on  ADPCM  [3.4]  but  none  on  Its  use  with  TASI. 
Operating  with  only  speech  signals,  computer  simulation  has  shown 
that  ADPCM/TASI  can  achieve  speech  compression  ratios  of  up  to 
3:1  before  any  degradation  is  perceived  by  a normal  user. 

4.1  SPEC 

The  SPEC  (Speech  Predictive  Encoding  Comjnunlcat Ion)  system 
was  conceived  In  1972  by  COMSAT  to  im.prove  the  efficiency  of 
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satellite  links.  It  was  originally  designed  to  accept  a I4.096 
Mb/s  TDM  frame  consisting  of  61t  PMC  voice  samples  every  125  usee, 
a European  standard.  The  output  rate  of  the  system  Is  2.0*18  Mb/s 
yielding  a bit  rate  compression  ratio  of  2:1.  Subsequent  proto- 
type models  of  the  system  could  achieve  a 2:1  compression  with 
any  number  of  Input  PCM  voice  channels,  e.g..  It  could  compress 
U8  channels  down  Into  the  bit  rate  of  21».  The  SPEC  system  Is 
described  In  [1]  and  field  test  results  over  INTELSAT  satellite 
links  are  contained  In  [2].  Other  field  tests  are  currently  In 
progress.  In  what  follows  we  suimiiarlze  the  configuration  of  the 
SPEC  system  when  used  to  compress  the  output  of  two  T1  terminals 
down  Into  a 1.5****  Mb/s  digital  signal  normally  required  by  one 
T1  terminal. 

A diagram  of  the  SPEC  transmitter  when  used  In  the  DDR 

application  Is  shown  In  Fig.  3.  The  two  T1  carrier  systems  shown 

contain  *18  voice  band  channels.  Two  of  these  channels  are 

permanently  assigned  to  data  and  are  not  affected  by  the  SPEC 

system*.  A time-shared  digital  voice  switch  continuously 

determines  which  of  the  **6  Input  voice  channels  contain  active 

speech  In  order  to  prevent  unnecessary  transmission  of  noise 

during  silent  Intervals.  The  zero-order  predictor  examines  each 

sample  from  each  Input  trunk  and  sends  only  unpredictable  PCM 

samples  along  with  a sample  assignment  word  to  the  output 

multiplexer.  A simple  parity  check  is  generated  by  the  error 

control  unit  to  protect  the  assignment  word.  The  asslgn-ment 

recirculation  unit  rotates  the  priority  of  service  to  the 

•Actually  for  the  DDR  application  from  2 to  6 channels  are 
dedicated  to  data.  This  explanation  assumes  only  two  data 
channels,  however. 
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Input  trunks.  Frame  and  assignment  sequence  synchronization 
codes  are  also  multiplexed  Into  an  output  TDM  frame.  The 
receiver,  shown  In  Fig.  U,  demultiplexes  the  received  frame  and 
performs  frame  and  sequence  synchronization.  Following  an  error 
check,  the  sample  assignment  word  Is  used  to  direct  unpredictable 
PCM  samples  to  the  correct  output  circuits.  Predictable  samples 
are  regenerated  by  a zero-order  predictor  located  In  the 
reconstruction  decoder. 

The  operation  of  the  SPEC  encoder  Is  Illustrated  In  Fig.  5. 

The  predictive  frame  memory  (Ppri)  stores  ^6  8-blt  words  (Pj^,P2, 

•••,  Pijg),  which  are  used  as  the  prediction  values  for  the 
corresponding  current  PCM  samples  (Sj^,S2, ‘ • ,S^g)  originating 
from  the  two  T1  carrier  lines.  The  current  samples  from  each 
trunk  Sj^,S2,  • • • along  with  the  corresponding  prediction  values 

Pi , P2 , • • • , Pjjg  processed  sequentially  by  the  predictor /detector. 

For  the  l'"^  trunk,  the  predictor/detector  compares  the  present 
PCM  samples  to  Its  corresponding  prediction  value  P^^.  If  the 
difference  between  and  P^^  Is  less  than  or  equal  to  a fixed 
threshold  6 then  Sj^  Is  not  transmitted  and  the  value  Pj^  remains 
In  the  1*^  location  of  the  PFM.  At  the  same  time,  a logic  0 Is 
written  Into  the  1^^  bit  position  of  a it6-blt  word  called  the 
sample  assignment  word  (SAW),  which  Is  part  of  the  transmitted 
TDM  frame.  Conversely,  If  the  difference  between  Sj  and  P^ 
exceeds  S,  then  Sj^  must  be  transmitted  and  therefore  Is  transferred 
to  one  of  the  8-blt  time  slots  Tj  of  the  output  frame.  The  value 
of  Sj  also  replaces  the  value  of  In  the  1^^  location  of  the 
PFM  and  a logic  1 is  written  In  . 
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The  output  TLiM  frame  shown  In  Fig.  5 consists  of  two  8-blt 
time  slots  reserved  for  data,  16  8-blt  time  slots  reserved  for 
transmission  of  the  unpredictable  voice  samples  and  eight  8-blt 
time  slots  that  contain  the  l<6-fclt  SAW.  The  output  frame 
Is  therefore  composed  of  2*)  8-blt  time  slots,  transmitted  In 
125  us,  resulting  In  an  output  rate  of  I.5UH  Mb/s.  This  repre- 
sents a bit-rate  reduction  of  exactly  2:1.  During  each  125  ys 
Interval,  a new  TDM  Input  frame  Is  delivered  to  the  SPEC  system 
and  a TDM  output  frame  Is  constructed  and  sent  to  the  receiver. 
Therefore,  during  any  125  us  Interval,  as  many  as  16  of  the  1(6 
voice  trunks  can  be  updated  at  the  receiver.  The  SAW  contains 
up  to  16  logic  I's,  which  are  distributed  over  1(6  bits  to  direct  the 
unpredictable  samples  that  are  transmitted  In  the  time  slots  Tj 
to  the  correct  outgoing  trunks  at  the  receiver. 

This  operation  of  the  SPEC  receiver  Is  illustrated  In  Fig.  6* 

As  the  transmit  frame  enters  the  receiver,  the  1(6  bit  SAW  is 
delivered  to  a sample  assignment  distributor,  which  directs  the 
new  samples  contained  In  time  slots  iTj . • • • ,T^g  to  the  correct 
location  In  the  FFM.  In  this  manner,  up  to  16  of  the  samples 
stored  In  the  receiver  PFM  are  updated  every  125  us.  All  samples 
not  updated  remain  the  same.  All  of  the  samples  contained  In 
the  PFM  are  then  passed  on  to  the  correct  channel  decoders  (In 
the  1,2, •••,^6  sequence)  every  125  us  to  reconstruct  the  speech 
signals . 

SPEC  Computer  Simulations 

In  order  to  access  the  performance  of  the  SPEC  system  In 
the  DDF  application  we  wrote  computer  programs  to  simulate  Its 
performance.  These  programs  were  run  on  an  ADAGE  AGT/30  computer 
and  objective  measurem.ents  made  of  SPEC  performance.  Audio  tares 
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were  also  made  to  demonstrate  system  performance.  These  tapes 
form  a part  of  this  report.  The  program  used  Is  called  SPEC3. 
A listing  of  this  program  Is  contained  In  Appendix  A of  this 
report. 

We  simulated  three  SPEC  systems  for  the  DDR  application. 
Their  properties  are  summarized  In  Table  7 below. 


SPEC 

System 

Active  Speech 
Channels 

Permanent 
Data  Channels 

Time  shared 
8-blt  speech 
slots/frame 

4 6/2 

46 

2 

16 

42/6 

42 

6 

12 

16/0 

16 

0 

16 

Table  7.  SPEC  Systems  Simulated 

The  *l6/2  system  has  46  active  (off-hook)  speech  channels  and 
two  data  channels.  "The  46  speech  channels  share  16  8-bVt  trans- 
mission slots.  Similarly  the  42/2  system  has  42  active  speech 
channels,  six  data  channels  and  12  8-blt  transmission  slots. 

The  l6/0  system  has  only  16  active  talkers  all  of  which  are 
assigned  to  an  8-blt  transmission  slot.  The  16/0  system  shows 
the  maximum  performance  of  SPEC  during  low  traffic  conditions. 

Its  performance  is  Identical  to  that  of  S-tlt  ij”255  PCK.  We  did 
not  consider  signaling,  supervision  or  alarm  arrangements  for 
these  systems.  The  46/2  system  Is  the  one  described  In  the 
previous  section.  The  SPEC  3 program  simulated  systems  which 
agreed  In  every  respect  except  two  with  the  system  described  in 
[1].  These  two  exceptions  are  outlined  in  Appendix  B.  Both 
changes  Improved  the  performance  of  the  system. 

The  SPEC  3 program  contains  subroutines  which  simulate 
traffic  loading  statistics. 
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The  following  eleven  soeech  segments  v;ere  used  as  source 
material  for  the  simulations: 

1.  Nine  rows  of  soldiers  stood  In  line,  (male  #1) 

2.  The  beach  Is  dry  and  shallow  at  low  tide.  (male 

3.  The  pipe  began  to  rust  while  new.  (female  speaker) 

I*.  Open  the  crate  but  don’t  break  the  glass.  (female  speaker) 

5.  Act  on  these  orders  with  great  speed.  (female  speaker) 

6.  Add  the  sum  to  the  product  of  these  three,  (male  #2) 

7.  Thieves  who  rob  friends  deserve  Jail.  (male  #2) 

8.  The  hog  crawled  under  the  high  fence.  (male  #2) 

9.  Cats  and  dogs  each  hate  the  other.  (male  1*3) 

10.  Open  the  crate  but  don't  break  the  glass,  (male  1I13) 

11.  Thieves  who  rob  friends  deserve  Jail.  (male  #3) 

Segments  1 and  2 were  made  using  a high  quality  microphone. 
Segments  3-11  are  copies  of  sentences  recorded  at  Bell  Telephone 
Laboratories  through  a 500  type  telephone  set.  These  sentences 
were  filtered  to  2^00  Hz  and  sampled  at  8 KHz.  Each  sample  value 
was  converted  Into  a 15-blt  number  using  a uniform  quantizer. 
These  15-blt  words  were  stored  on  disk  In  the  AGT/30  computer 
and  used  as  source  material  for  the  encoding  systems  simulated 

on  this  contract.  Table  8 shows  some  of  the  details  of  these 

sentences  when  run  through  an  8-blt  u-255  companded  PCM  system. 

Since  the  sampling  rate  Is  8 KHz  the  number  of  samples  shows 

that  each  segment  Is  from  two  to  three  seconds  long.  The  slgnal- 

to-quart Izlng-nolse  ratios  of  about  37dB  verify  that  the  8-blt 

U“255  PCM  system  Is  Ideal.  The  S/N*  ratios  are  almost  Identical 

for  all  speech  segments.  The  signal  power  merely  reflects  the 

•In  this  report  the  term  S/N  refers  to  the  ratio  of  signal  power 
to  unweighted  noise  power.  In  general  "C"  message  weighting  of 
the  noise  adds  2 . 9dB  to  this  ratio. 
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Table  8.  Source  r^terlal  for  computer  simulations 


Speech 

Segment 

No,  of 
Samples 

Slg  Power 
In  dBm 

Power  relative  to  a 
full  load  sinusoid 

1 

23219 

37.2 

-16.8 

-13.9dB 

2 

23057 

37. It 

-16,3 

-13.3dB 

3 

15765 

37. >> 

-16.9 

-13.9dE 

It 

19025 

37.2 

-15.5 

-12.5dB 

5 

17609 

37.2 

-17.1 

-lit.ldB 

6 

21^165 

37.3 

-19.7 

-l6.7dE 

7 

21151 

37.3 

-19.7 

-l6.7dE 

8 

20629 

37.5 

-18.3 

-15.3dB 

9 

16779 

37.2 

-20.8 

-17.BdB 

10 

16931 

37.5 

-20.0 

-17.0dB 

11 

17361 

37.2 

-19.5 

-l6.5dE 

computer  scaling  of  the  sample  values.  The  column  entitled 
"Power  relative  to  a full  load  sinusoid”  shows  how  the  stored 
speech  segment  Is  matched  to  the  simulated  PCM  system.  The 
power  In  the  signal  can  be  changed  programmatically  before  passing 
the  signals  through  any  of  the  computer  simulations. 

Figure  7 shows  how  the  performance  of  the  three  SPEC  systems 
varies  with  Input  power  level.  The  curves  show  the  dynamic  range 
of  about  35dB  typical  of  u«255  PC.V  systems.  The  l6/0  system  Is 
operating  like  an  ideal  PCM  system  while  the  U6/2  and  >*2/6 
systems  have  naxlr.uir.  S/N  ratios  of  about  29  and  23d3  respectively. 

Figure  8 shows  how  the  performance  of  the  SPEC  systems 
changes  with  the  number  of  off-hook  channels.  The  *16/2  system  has 
16  8-blt  transmission  slots  per  frame  while  the  ^2/6  system  has 
only  12.  This  Is  the  only  figure  In  which  the  systems  are  not 
fully  loaded.  The  fact  that  these  computer  simulations  don't  lie 
on  a smooth  curve  and  seem  slightly  erratic  Is  due  to  the  fact 
that  each  point  represents  a simulation  of  only  5 seconds  of  speech. 
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Figure  7.  S/S  vs.  signal  level  for  th.^ee  SPEC  systems  usl.nr 
speech  segments  3 and  7. 
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S/N  ratio  vs.  number  of  off-hook  channels  for  the 
1)6/2  and  1)2/6  SPEC  systems.  Input  Is  speech 
segments  3 and  7. 


For  the  M6/2  system  once  the  loading  goes  down  to  about  3**  off-hook 
channels  the  system  performance  approaches  that  of  8-blt  PCM. 

For  the  1(2/6  system  the  loading  must  be  about  28  off-hook 
channels  before  the  performance  approaches  that  of  8-blt  PCM. 

Figure  9 shows  the  S/N  vs.  probability  of  error  for  the  three 
SPEC  systems. 

U.2  ADPCM/TASI 

Five  bit  DPCM  with  an  adaptive  quantizer  (ADPCM)  achieves  a 
S/N  ratio  for  speech  cf  about  31<1B  [ll].  For  8-blt  PCM  with  a 
D-channel  bank  the  S/N  ratio  Is  about,  37.5dB.  In  an  operational 
environment,  however,  PCM  S/N  ratios  of  33-36tlB  are  the  rule. 
Cumr.eskey,  ^.al^.  [3]  have  shown  that  l(-blt  ADPCM  Is  (almost) 
subjectively  equivalent  to  7-blt  PCM  (u”100).  The  subjective 
performance  of  ADPCM  Is  better  than  Its  S/N  ratio  indicates. 
Cumraeskey's  experiments  and  others  show  that  ADPCM  can  transmit 
speech  with  2 to  3 fewer  bits  per  sample  than  PCM.  In  other 
words,  a 5-  or  6-blt  ADPCM  should  give  the  same  subjective  per- 
formance as  an  8-blt  u*255  PCM  system.  A 5-blt  ADPCM  system 
would  probably  be  slightly  worse  than  8-blt  PCM  and  a 6-blt  DPCM 
system  would  be  slightly  better.  These  conclusions  have  been 
verified  by  computer  simulations  run  here  at  NCSU  using  the 
computer  program  JADPCM. 

When  used  with  the  TASI  system  described  In  Chapter  III  of 
this  report  It  Is  necessary  for  the  ADPC.M  system  to  be  able  to 
drop  bits  during  periods  of  high  activity.  Thus  the  ADPCM/TASI 
system  may  have  word  lengths  of  8,  7,  6,  5 or  even  fewer  bits 
depending  on  the  number  of  active  talkers.  As  with  SPEC  a voice 
switch  determines  whether  or  not  speech  Is  present  on  each 
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Input  channel.  If  It  Is,  the  ADPCM  output  Is  connected  to  a 
transmission  channel.  If  not,  no  sample  Is  transmitted.  For  each 
frame  the  system  must  determine  how  many  bits  are  available  for 
each  sample  value.  For  example,  for  the  U2/6  system  there  are 
132  bits  per  125  usee  frame.  If  there  are  2^1  active  channels 

then  each  channel  Is  entitled  to  either  5 or  6 bits  because 
1?2  1 

» Sj-  . Thus  12  channels  will  be  assigned  L»5,  and  the  other 
12  will  have  L"6.  The  program  JADPCM  simulates  the  ADPCH  system 
described  In  [3]  along  with  the  TASI  system  described  here.  The 
loading  statistics  were  simulated  programmatically  as  In  the  SPEC3 
program.  Further  details  of  the  operation  of  ADPCM/TASI  are 
covered  In  the  following  chapter  of  this  report. 

Computer  simulation  of  the  ADPCW/TASI  system 

The  ADPCM/TASI  system  has  been  simulated  using  the  computer 
program  HCDAQ2.  Figure  10  shows  the  dynamic  range  of  the  system. 

It  Is  somewhat  broader  than  that  of  standard  PCM.  The  adaptive 
algorithm  used  was  that  described  In  [3]-  The  dynamic  range  of  the 
system  Is  a function  of  that  algorithm.  This  adaptive  algorithm 
also  limited  the  maximum  S/N  ratio  to  about  32dB.  This  could  be 
Increased  to  37dB. 

Figure  11  shows  the  performance  of  the  system  vs.  number  of 
off-hook  channels  for  speech  segments  3 and  10.  The  16/0  system 
has  the  same  S/N  ratio  as  the  46/2  system  and  Is  not  plotted. 

Figure  12  shows  how  the  expected  performance  of  the  system 
varies  as  a function  of  the  probability  of  error  on  the  line. 

The  MCDAQ2  program  does  not,  at  this  time,  have  the  capability  of 
Introducing  errors.  We  do  have  a program  JADPCM  which  has  the 
same  ADPCM  algorithm  wlt.hout  the  variable  step  size  capability 
and  JADPCM  can  simulate  errors  on  the  line.  In  Pig.  12  we  show 
the- performance  of  ADPCM  with  4 and  5 bits  per  sample  value. 
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Figure  10.  S/N  vs.  signal  level  Tor  ADPC^/TASI.  Points  show 
results  from  program  MCDAQ  using  speech  segments 
3 and  10.  This  figure  applies  to  the  16/0,  <t6/2 
and  t2/£  systems. 
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Figure  11.  ADPCM/TASI  performance.  S/N  vs.  off-hook  Input 

channels  for  the  *16/2  and  ll2/6  systems  using  speech 
segments  3 and  10. 
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Igure  12.  Probability  of  error  vs.  S/N  ratio  for  ADPCW.  Points 
show  results  from  computer  simulation  uslnn  orogram 
JADPCM.  Also  shown  are  results  projected  for  the 
ADPCK/TASI  system.  The  '-*6/2,  1*2/6  and  16/0  systems 
would  have  Identical  performance. 


From  these  results  we  have  extrapolated  the  expected  error 
performance  of  the  ADPCM/TASI  system.  There  Is  no  essential 
difference  In  performance  between  the  46/2,  42/6  and  l6/0  systems. 
The  S/N  plotted  In  Fig.  12  shows  that  the  ADPCM/TASI  system  Is 
more  sensitive  to  errors  In  the  transmission  medium  than  SPEC  or 
standard  PCM.  This  Is  due  to  the  adaptive  feature  of  the  ADPCM 
system.  When  an  error  causes  the  decoder  step  size  to  be  dif- 
ferent from  that  at  the  encoder,  a considerable  discrepancy 
between  the  Input  and  output  signals  results.  This  discrepancy 
Is  a difference  In  signal  level  and  Is  not  nearly  so  degrading  as 
noise  of  equal  power.  Thus  the  subjective  performance  of  the 
system  Is  considerably  better  than  the  S/N  ratio  In  Fig.  12 
would  Indicate.  The  degradation  exhibits  Itself  subjectively  as 
a random  change  In  the  speech  level. 

4.3  NIC/TASI 

Nearly  Instantaneous  Companding  (NIC)  Is  a scheme  for 
processing  15-segment  u-law  or  IS-segment  A-law  PCM  to  achieve 
a reduced  bit  rate  [8].  This  scheme  Is  a type  of  APCM.  Speci- 
fically, the  NIC  coder  (see  Fig.  13)  groups  the  n-blt  u-law  PCM 
encoded  words  at  its  Input  Into  groups  of  N and  looks  for  the 
sample  havl.ng  the  largest  magnitude  In  the  group.  It  then  re- 
encodes each  sample  In  the  group  to  approximately  n-2  bit  uniform 
quantization  with  overload  at  the  top  of  the  chord  of  the  maxlmuni 
magnitude  sample.  These  n-2  bit  uniform  words  are  transmitted 
to  a complimentary  SIC  decoder  that  regenerates  n-blt  codes  In 
the  original  companding  law.  Since  three  scaling  bits  giving 
the  chord  of  the  maximum-magnitude  sample  must  also  be  transmitted 
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Figure  13.  Processing  8-blt  PCM  with  the  NIC  algcrltnn. 
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to  the  decoder,  the  resulting  bit  rate  Is  f.(n-2+3/N)  where  f 
Is  the  sampling  rate  of  8 KHz.  This  bit  rate  Is  comparable  to 
that  of  n-2  bit  PCM  for  even  modestly  large  N.  For  the  DDR 
system  we  use  ;,'-10.  The  dynamic  ra.nge  of  the  original  PCM 
companding  law  Is  retained  through  the  ad  jtlve  overload  point 
as  Is  the  relative  PCM  Insensitivity  to  signal  statistics. 

In  the  NIC/TASI  system  the  193-blb  frame  would  be  divided 
as  follows:  one  bit  for  fra.mlng,  12  bits  for  voice  switch  and 
signaling  Information,  180  bits  would  be  shared  amor.g  the  dedi- 
cated data  channels  v2  to  6 8-blt  words  per  frame)  and  the  speech 
channels  (1(6  to  112).  The  1*6/2  system  would  have  16  bits  per 
frame  dedicated  to  data  and  161*  tits  would  be  shared  among  the 
1*6  speech  channels.  The  1*2/6  system  would  have  1*8  blts/frane 
dedicated  to  data  and  132  bits  shared  among  the  i*2  speech 
channels.  This  Is  the  same  framing  arrangement  shown  In  Fig.  2. 

For  an  NIC  system  with  N”10  the  number  of  bits  per  speech 
sample  Is,  on  the  average,  6+^  • 6.3.  ?cr  the  i*6/2  system  the 
number  of  bits  available  per  frame  for  speech  samples  Is  16**. 

If  26  channels  are  active  the  number  of  bits  required  Is 
26  X 6.3  - 163.8  — thus  If  26  or  fewer  channels  are  active  the 
NIC/TASI  system  will  operate  as  an  NIC  system.  If  more  than  26 
channels  are  active  then  some  bits  must  be  dropped  frcm  the  6-blt 
NIC  word.  The  probability  of  this  happening  Is  the  probability 
that  27  or  more  channels  are  active  at  the  same  time.  The 
probability  of  this  happening  Is  E(27,<*6, . 33)  ■ .003-  Thus  for 
the  1*6/2  system  a bit  would  seldom  be  dropped  and  the  NIC/TASI 
system  should  perform  like  standard  NIC.  For  the  U2/b  system 
there  are  132  bits  per  frame  and,  since  6.3  x 20  • 126,  21  channels 
would  have  to  be  active  before  bit  dropping  Is  necessary. 
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The  probability  that  a bit  v/ould  need  to  be  dropped  on  a given 
speech  sample  Is  B(21,U2,.33)  - .038.  Since  this  Is  such  a 
small  percentage  of  the  time  It  would  have  little  effect  on  the 
S/N  ratio  and  again  the  performance  of  the  system  would  be  very 
nearly  like  that  of  NIC  without  companding. 

Unfortunately  there  Is  very  little  published  Information 
about  NIC.  We  became  aware  that  NIC  Is  a possible  candidate 
system  for  the  DCR  In  September  when  the  August  1976  Issue  of  the 
IEEE  Transact  Ions  on  Comnuncat Ions  became  available  [8].  No 
subjective  tests  on  the  system  have  been  published  and  the 
performance  of  the  system  In  the  presence  of  errors  Is  not 
documented.  Thus  at  the  present  tlm.e  [11/29/76]  we  know  only 
that  It  may  be  a promising  system  for  bit  rate  reduction.* 

it  .It  Comparison  of  Candidate  Systems 

In  this  section  we  discuss  the  relative  m.erlts  of  the  three 
candidate  systems,  SPEC,  NIC/TASI  and  ADPCM/TASI  and  describe 
the  reasons  for  choosing  the  ADPCK/TASl  system.  In  preference  to 
the  other  two. 

We  have  compared  these  systems  based  on  the  five  criteria 
listed  In  Table  9-  In  each  colunr.  of  the  table  the  performance 
of  each  system  Is  rated  relative  to  the  other  two  — a rating  of 
1 means  most  desirable  and  3 Is  least  desirable. 

Table  9.  Comparison  of  candidate  TDR  systems 

Attribute  ADPCH/TASl  NIC/TASI  SPEC 

Equipment  complexity  3 I I 

Performance  with  light  traffic  load  1 11 

Performance  with  heavy  traffic  load  1 32 

Performance  In  presence  of  errors  2 11 

Possibilities  for  further  corasresslor  1 2 2 

•Since  the  text  of  this  report  was  written  we  have  simulated  the  NIC 
system  and  found  Its  performance  to  be  less  than  that  of  SPEC  for 
the  DDR  application. 
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On  the  basis  of  eo.ulpment  complexity  the  ADPCM/TASI  system 
Is  more  complex  than  the  other  two  and  Is  rated  2 while  SPEC 
and  .VIC/TASI  seem  to  be  equally  complex  and  are  rated  1.  Each 
system  requires  essentially  the  same  voice  switch  hardware  and 
this  will  be  a major  cost  of  the  system.  The  ADPCV  system  must 
contain  hardware  for  the  differential  operation,  which  Is 
difficult  with  a companded  Input  and  step  size  adaptation.  This 
Is  t.he  primary  reason  It  Is  more  complex  than  the  other  two 
systems.  The  NIC/TASI  system  must  contain  a procedure  for 
selecting  the  largest  sample  In  a group  of  ten  and  some  sort  of 
buffering  and  multiplexing  arrangement  to  transmit  the  3-blt 
chord  Information  every  10  words.  The  SPEC  system  must  contain 
a priority  selection  procedure  and  hardware  for  the  variable 
threshold  arrangement.  All  three  systems  must  extract  signaling 
from  the  PCM  words  every  sixth  frame  and  ree.ncode  these  bits 
Into  a signaling  channel  which  must  be  multlclexed  into  the  bit 
stream.  All  three  systems  must  derive,  multiplex  and  transmit 
channel  assignment  information.  In  addition  to  this  the  .NIC/TASI 
and  ADPCM/TASI  systems  must  protect  their  channel  assignment 
Information  with  an  error  correcting  code.  The  alarm,  fram.l.ng 
and  interfacing  hardware  are  alike  for  all  three  systems.  If 
such  systems  were  ultimately  realized  through  1.SI,  the  cost  of 
all  three  systems  would  be  essentially  the  same  even  if  there 
are  differences  In  their  complexity. 

The  performance  of  the  systems  during  heavy  load  could  be 
quite  different.  Our  computer  simulations  show  that  ADPCM/TAEI 
has  no  perceivable  degradation  even  with  all  channels  off-hook 
for  both  the  <*6/2  and  the  <<2/6  versions  of  the  DDP.  The  SPEC 
system  simulations  show  that  the  noise  Is  noticeable  with  the 


^6/2  system  and  objectionable  with  the  1*2/6  system.  Objective 
S/N  comparisons  shown  In  Fig.  7 for  SPEC  and  Fig.  10  for  ADPCM/ 
TASI  verify  the  performance  advantage  of  ADPCM/TASI. 

The  NIC/TASI  system  has  recently  been  simulated  here  at 
NCSU  and  at  Bell  Telephone  Laboratories  [8j.  Their  con- 
clusion Is  that  the  objective  S/N  ratio  of  NIC  Is  about  3dB 
worse  than  could  be  provided  by  ADPCM.  Other  simulations  of 
ADPCM  [**]  v/hlch  have  bee.n  verified  here  show  that  ADPCM  Is 
subjectively  about  3ilB  better  than  Its  objectively  measured  £/N 
ratio  Indicates.  In  comparison  to  PCM  for  example  an  ADPCM 
system  with  a S/N  ratio  about  3<1B  below  that  of  a PCM  system  will 
be  subjectively  equivalent  to  the  PCM.  system..  E.xtrapolatlon  from, 
these  subjective  observations  shows  that  ADPCM/TASI  will  sound 
6dB  better  than  NIC/TASI  even  though  Its  S/N  ratio  may  only  be 
3dB  higher.  This  Is  significant  for  a heavily  loaded  system. 

For  this  reaso.n  we  rated  ADPCM/TASI  "1".  NIC/TASI  "3"  (although 
little  Is  known  about  It)  and  SPEC  ”2". 

The  performance  of  the  system.s  during  a light  load  are 
subjectively  equivalent  - noise  will  not  be  perceived  on  any 
of  them. 

Or.  performance  In  the  presence  of  errors  we  rated  SPEC  "1", 
ADPCM/TASI  "2"  and  NIC/TASI  "1".  Figs.  9 and  12  show  that  the 
objectively  measured  S/N  ratio  for  ADPCM/TASI  Is  slightly  worse 
than  the  performance  of  SPEC.  The  error  performa.nce  curves  for 
NIC  have  not  been  published.  It  has  been  surmised,  however, 
t.hat  NIC  has  an  error  performance  like  that  of  PCM  [3].  If  this 
Is  true  then  NIC/TASI  may  have  better  performance  In  the  presence 
of  errors  than  either  of  phe  other  two  competing  systems. 
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On  possibilities  for  further  coir.presslon  ADPCK/TASI  Is 
clearly  better  than  the  other  two  systems.  Computer  simulations 
of  this  system  were  made  for  a compression  ratio  of  3:1  and  it:l. 
The  3:1  compression  system  compressed  three  T1  systems  down  Into 
the  bit  stream  normally  occupied  by  one  system.  No  degradation 
Is  perceived.  The  £/H  ratio  Is  greater  than  30dB.  Simulations 
were  also  run  for  : 1 compression  — some  noise  can  be  detected 
but  It  Is  not  objectionable.  Cubjectlve  comparisons  with  SPEC 
show  that  the  ADPCK/TASZ  Is  clearly  superior.  The  NIC/TASI 
system  was  rated  "3"  In  this  category.  Even  though  very  little 
Is  known  about  this  system  we  feel  It  would  perform  better  than 
SPEC  because  SPEC  would  tend  to  freeze  out  with  very  heavy  loading 
while  NIC/TASI  will  always  transmit  some  bits  for  each  channel 
during  every  sampling  Interval. 

Making  a choice  between  these  three  candidate  systems  Is 
not  easy.  They  each  have  favorable  aspects.  However, 
because  of  Its  good  performance  under  heavy  load  and  Its  potential 
for  even  further  compression  we  feel  that  the  ADPCM/TASI 
system  Is  the  test  choice  for  the  CDR  application. 
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5.  PRELIMINARY  DESIGN  OF  AN  ADPCM/TASI  SYSTEM 


5 . 1 Block  Dlajrrams 

Figures  1,  2 and  3 Illustrate  a preliminary  design  of  an 
ADPCM/TASI  system  for  the  DDR  application.  For  the  purposes  of 
this  discussion  we  assume  that  the  system  Is  designed  to  transmit 
up  to  six  6^  K blt/sec  data  channels. 

Figure  2 shows  the  frame  transmitted  by  the  ADPCM/TASI 
system.  As  in  PCM  each  frame  consists  of  193  bits  which  occur 
every  125usec.  This  125usec  Is  the  Inverse  of  the  PCM  sampling 
rate  of  8 KHz.  One  bit  every  frame  Is  devoted  to  framing.  The 
remaining  192  bits  are  divided  Into  three  separate  channels  — 
a signaling  channel  of  8 bits,  a channel  assignment  signal  of  2 bits, 
data  channels  consisting  of  six  8-blt  words  for  a total  cf  i48  bits 
and  speech  transmission  bits  consisting  of  the  remaining  132  bits. 

For  the  purposes  of  this  discussion  we  have  assumed  that  there 
are  six  8-blt  data  channels  per  frame  as  shown  In  Fig.  Kb). 

Actually  this  Is  adjustable  from  0 to  six  8-blt  words  per  frame. 

The  signaling  channel  and  channel  assignment  signals  always 
contain  8 and  4 bits  per  frame,  respectively.  The  number  of 
speech  transmission  bits  can  be  as  high  as  180  when  there  are  no 
dedicated  data  channels  or  as  low  as  132  when  there  are  six 
data  channels. 

The  encoder  of  Fig.  14  shows  that  the  two  T1  Inputs  have  a 
common  clock.  If  this  common  clock  is  not  provided  then  the  T1 
Input  lines  rust  be  buffered  at  the  Input  of  the  encoder  and 
bit  stuffing  techniques  must  be  used  to  Insure  that  the  1.5*^**  m3 
lines  are  In  sync,  '’bis  can  be  done  tut  hopefully  Is  not 
necessary. 
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The  "signaling  extraction"  box  renoves  the  signaling  In- 
formation from,  all  speech  channels  and  forns  a separate  S-blt 
per  fram.e  signaling  channel.  The  "tlm.lng  extraction"  function 
deter-.lnes  frame  and  bit  synchron 1 cat  Ion . Although  not  shown 
on  the  diagram,  these  tiring  signals  are  passed  to  the  ot.ner 
parts  of  the  system  which  need  then. 

The  "external  priority  selection"  function  Is  a set  of 
switches  or  signals  which  determ.lne  which  of  the  A8  Ir.out  speech 
channels  should  bypass  the  ADPC.V  operation.  The  data  channels 
would  be  sc  designated  as  could  any  other  ch.ar,nel  for  which  the 
ADPCy  operation  night  prove  unsat Isfactcry . The  "external 
priority  selection"  controls  switch  S.  which  feeds  the  data 
signals  to  a data  buffer  bypassing  the  ADPCM  ope.-atlon. 

Essential  to  the  operation  of  any  TASI-type  system,  are  the 
voice  switches,  shown  In  Plg.lA,  that  car  recognize  the  presence 
or  absence  of  speech  In  a noisy  environment.  This  task  Is  made 
quite  difficult  by  the  fact  that  many  speech  sounds  are  like 
noise.  All  voice  switches  Introduce  a certain  amount  of  clipping 
at  the  beginning  of  a speech  burst.  Speech  clips  longer  than 
about  50  m sec  cause  objectionable  mutilation  of  Initial  frica- 
tive, nasal  and  plosive  consonants  [li!].  The  state-of-the-art 
In  voice  switch  design  Is  currently  represented  by  the  voice 
switch  designed  for  the  SPEC  system  by  CCMSAT.  It  operates  off 
of  an  8-blt  u*255  PCM  bit  stream.  The  starting  point  for  our 
voice  switch  will  be  this  COKSAT  design  which  Is  described  In 
[17].  The  number  of  voice  switches  on  Is  passed  to  "co.r.pute  L" 
which  computes  the  nu.mber  of  bits  L available  for  each  speech 
channel. 
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The  box  marked  "control"  has  as  Its  Input  the  number  L,  the 
voice  switch  Information  for  each  speech  channel,  and  the  priority 
Information.  From  these  Inputs  It  creates  the  control  signal  for 
switch  S^.  This  switch  selects  the  channel  whose  voice  switches 
are  on  and  assigns  them  a transmission  slot  L bits  wide. 

This  control  Information  contains  the  channel  assignment.  It  Is 
also  encoded  using  a rate  j error  correction  code  and  multiplexed 
onto  the  l.b***!  Fb/btransnltted  bit  stream.  The  error  correction 
code  Is  needed  to  protect  this  control  Information  against 
errors  In  the  transmission  medium. 

The  "8-blt  log  to  13-blt  linear  converter"  converts  the 
8-tlt  u*255  PCM  Input  Into  a 13-blt  linear  code.  This  Is  done  so 
that  the  AtPCM  system  (enclosed  by  a dashed  box  In  Fig.  !<))  can 
operate  with  a linear  digital  code.  The  smallest  step  In  an  8-blt 
'.j*255  companded  code  Is  equivalent  to  the  smallest  step  In  a 
13-blt  linear  code.  This  code  conversion  Is  relatively  straight- 
forward [l9]. 

The  quantizer  In  the  ADPCM  system  Is  the  "13  bit  to  L bit 
converter".  The  feedback  loop  consists  of  the  "L  bit  to  13-tit 
converter",  the*13-blt,  **8  word  menory'and  the  multiplier  M. 

If  we  choose  M“.S75,  a reasonable  value,  this  r.ultlpller  becomes 
a slr.ple  shift  and  add.  The  "step  size  adaptation"  procedure  Is 
tr.at  reccr.T.ended  by  Jayant  [3]- 

The  output  of  the  ADPC.M  system  Is  stored  In  the  "8-blt  *»8 
word  memory  for  speech".  The  I,  bit  speech  samples  are  held  there 
until  the  switch  connects  them  to  an  appropriate  transmission 
channel.  The  data  sar.ples  are  held  In  the  "8-blt  6 word  memory 
for  data"  until  the  switch  S,  puts  them  Into  a data  slot. 
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Figure  15  Is  a block  diagram  of  the  ADPCM/TASI  decoder. 

Its  operation  Is  the  Inverse  of  the  encoder  and  need  not  be 
discussed  further. 

5.2  Hardware 

It  Is  quite  difficult  to  estimate  the  hardware  requirements 
from  the  block  diagrams  of  Figs.  llJ  and  15.  The  chip  count  Is 
dependent  on  the  detailed  design  which  has  not  yet  been  done. 

It  may  be  desirable  to  Impler.e.nt  the  encoder  and/or  decoder 
using  prograr.med  controllers.  The  scale  or  size  of  the  Inte- 
grated circuit  packages  to  be  used  has  not  yet  been  determined. 
Nevertheless,  we  estimate  that  the  chip  count  for  the  encoder  to 
be  approximately  150  to  300  chips.  For  the  decoder  we  expect  a 
chip  count  from  100  to  250.  These  estimates  include  the  T1 
hardware  to  be  purchased  for  the  "signaling  extraction",  "timing 
extraction"  and  "alarms"  fu.nctlons. 
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6.  AUDIO  TAPE  DEHCNSTRATIOH 

An  audio  tape  was  made  to  demonstrate  the  performance  of 
some  of  the  data  compression  systems  described  In  this  report. 
This  audio  tape  Is  net  furnished  with  this  report.  The  tape 
speed  Is  7-1/2  Inches/sec.  and  the  recording  Is  on  channel  1. 
Except  for  Parts  A and  J a duty  factor  of  .33  was  used. 

The  contents  of  the  tape  are  as  follows: 


.''art  A:  Part  A contains  the  11  speech  segments  used  as  source 

•material  for  the  computer  simulations.  These  sentences 
are  listed  In  section  IV-A  of  this  report. 

Part  E;  16/C  System  - no  errors 

Vale  speaker  - then  female 

1.  Original 
3.  ADPCM/TASI 

Part  C;  A6/2  System  - no  errors 

l^ale  sr -'aker  - then  female 

1.  Original 

2.  SPEC 

3.  ADPCK/TASI 


Part  D: 


“2/6  System  - no  errors 

Yale  speaker  - then  female 

1.  Original 

2.  SPEC 

3.  ADPCM/TASI 


Part  E:  Performance  of  SPEC  in  the  presence  of  errors  - l6/0 
system,  .'•'ale  speaker  only  but  with  1 repetition 


1.  Original  ^ 

2.  P(e;  - 10';; 

3.  P(e)  - 1C', 
A.  P(e)  - 10'^ 
5.  P(e)  - 10'*^ 
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Fart  F:  Perforrance  of  5-blt  JADFCM  In  the  presence  of  errors. 

This  Is  a lower  Unit  on  the  perfornance  of  ADPCK/TASI 


- 16/0  system.  Male  speaker  with  one  repetition. 

1.  Original  - 

2.  P(e)  - 10“? 

3.  P(e)  - 10', 
i*.  P{e)  - 10'^ 

5.  P(e)  - 10'“^ 

Part  G:  Performance  of  5-blt  JADPCM  compared  with  SPEC  with  a 

probability  of  error  of  P(e)  • 10  • Male  speaker  with 

one  repetition. 

1.  Original  - no  errors 

2.  SPEC  - P(e)  - lO-**  , 

3.  JADPCM  - P(e)  - 1C" 

Part  K:  Performance  of  5-blt  JADPCM  compared  with  SPEC  with  a 
F(e)  • lO"'.  Male  speaker  with  one  repetition. 

1.  Original  - no  errors 

2.  SPEC  - P(e)  - 10-3 

3.  JAD.^CM  - P(e)  « lO"-’ 

Part  I:  3:1  and  ^1:1  compression  with  ADPCfl/TASI: 

3 T1  systems  compressed  into  1.5^t  mB  and  T1  systems 
compressed  Into  1.5^**  mB  - .no  dedicated  data-  Male 
speaker  with  one  repet  It  lo.n. 

1.  Original 

2.  3:1  compression  (177-blt  frame) 

3.  I*:!  compression  (172-blt  frame) 

Part  J:  Co.T.parlsor.  of  ADPCM./TASI  with  SPEC  for  42/6  systems  with 
duty  factor  of  a«.45.  Male  speaker  with  one  repetition. 

1.  Original 

2.  SPEC 

3.  ADPCM/TASI 
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APPENDIX  A 


Computer  Programs 

A number  of  computer  programs  were  written  during  the 
course  of  this  study.  In  this  appendix  we  list  the  two  most 
Important  ones: 

1.  MCDAQ2  - Multiple  Channel  Differential  PCM  with  an 
Adaptive  Quantizer.  This  program  simulates  the  ADPCM/ 
TASI  system  described  in  this  report. 

2.  SPEC3  - Speech  Predictive  Encoding  Communication  system. 
This  program  simulates  the  SPEC  system  described  In 

this  report. 

These  programs  are  written  In  FORTRAN  IV  G Level,  Adage 
Version.  In  addition  to  these  programs  a large  number  of  support 
programs,  subroutines  and  Input/output  procedures  were  written. 

If  all  of  this  software  were  Included  here  It  would  double  the 
size  of  this  report.  All  of  the  software  Is  available  at  North 
Carolina  State  University. 
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PROGRAM  SPECS 

C SIMULATION  OF  COMSAT  SPEC  SYSTEM 
C WITH  ONE  CHANNEL  OF  REAL  SPEECH. 

C THE  REMAINING  CHANNELS  ARE  MODELED  BY  RANDOM 
C EXPONENTIALLY  DISTRIBUTED  VOICE  SWITCH  ON  AND  OFF 
C periods,  and  INDEPENDENT  SEQUENCES  OF  UPDATE  REQUESTS. 

C quantization  is  A-LAW,  13  SEGMENT,  8 BIT. 

C R.W.S.  2/25/76 
C VERS.  1,  REV.  C: 

C SOME  NEW  STATISTICAL  MEASUREMENTS  ADDED 

C AND  SNR  CALCULATIONS  MODIFIED  FOR  LESS  TRUNCATION  ERROR. 

C VERS.  1,  REV.  D,  3/2/76 

C FREEZEOUTS  & FREEZEOUT  ERROR  ACCUMULATED  AS 
C A FUNCTION  OF  ACTIVITY 

C VERS.  1,  REV.  E,  3/29/76 

C ADAPTIVE  THRESHOLD  ADDED  AS  AN  OPTION 


C VERS.  1,  REV.  F,  6/1/76... 

C ADAPTIVE  THRESHOLD  RULES  FOR  NO  = 16  OR  18 
C SOME  STATISTICS  REMOVED  FOR  SIMPLICITY 


C VERS.  1,  REV.  G,  6/1/76... 
C BIT  ERROR  SIMULATION  ADDED 


C VERS.  1,  REV.  H,  7/23/76... 

C ADAPTIVE  THRESHOLD  RULE  FOR  NO  = 12  ADDED. 


C SPEC3  IS  SPECl  MODIFIED  TO  USE  FUNCTIONS 
C PCMWD  &PCMDC  IN  PLACE  OF  Q255  & UQ255 
C JB0, 9/7/76 


INTEGER  RBUF{1664) ,WBUF(1664) ,MSLST(20) ,CNT(100) 
INTEGER  ACTIV(48) ,ADTHR(48) ,KZOA{48) 

INTEGER  DNT,DSGN 

INTEGER  PCMWD .DELTA , PCODE , DCODE , EXPRN 
FRACTION  PCMDC,SAMP,QSMP 

LOGICAL  END, STATE (100) ,OFFLG , PFLG , CH I F , ROTFG , ADTFG 
DATA  SPLNG/1 . 5/ , ALPHA/ .38/ ,QSCL/1 . / 

DATA  MSLST/20*0/ ,DELTA/0/,DT/125.E-6/ ,URNR/134217728./ 
DATA  PURDS/.87/ 
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C 
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OF  EXCEEDING 


ADAPTIVE  THRESHOLD  RULE  & PROBABILITIES 
THRESHOLD  VS.  THRESHOLD: 

INTEGER  ADTHR16(48) ,ADTHR18 (48) .ADTHR12(48) 

DIMENSION  PURDA16 (48) ,PURDA18 (48) ,PURDA12 (48) 

DATA  ADTHR12/12<^''  n 

DATA  ADTHR16/16’ 

DATA  ADTHR13/18’ 

DATA  PURDA12/12* .87 , . 76 . . 68 . 2* . 59 . 3* . 53 , 2 * . 48 , 27 * . 44/ 
DATA  PURDA16/16* .87 ,2* .76 ,.68, 4*. 59, 3*. 53, 22 *.48/ 

DATA  PURDA18/18* .87 ,2* .76,2* .68,4* .59,3*.53,19* .48/ 


C BIT  error  SIMULATION  VARIABLES: 

INTEGER  TCSB,PCMEP,SAW(48) ,ESAW(48) ,EBIT,ESEED 

LOGICAL  SBEFG,PBEFG,SERRF,PERRF 

LOGICAL  CRDFG 

ESEEU=1234567 

PE=. 000001 

ISEEU=7654321 

CRDFG=.F. 


760 

700 

730 


770 

/20 


750 


GO  TO  710 
WRITE (25.700) 

FORMAT ( 1 CARD  READ  OPTION') 

CALL  YESNO (CRDFG) 

1F{  .N.CRDF(3)  GO  TO  710 
ESEED=1234567 

READ(17.720}  NI  . NO . MSLST ( 1 ) , SPG, PE 

FORmAt(3I2,F10.5,f15.10) 

IF(PE.EQ.0)  PE=. 000001 
IF(NI.Eg.0i  GO  TO  740 
WRITE  (2  5/7  50)  NX  , NO,  MSLST  (1)  .SP(3,PE 
FgRMA-f  1iNp6t=  ,12,/,  O0TPCiT=  'i2, 
PROB  ERR=  ,F10.8) 


GO 


GAIN=  ■ ,F6.2,// 
TO  780 


/, 


SEGMENT= 


12,/, 


710 

900 

901 

902 
10 

903 

904 

913 

914 

915 

916 


WRITE (25,900) 
FORMAT ( 0NO.  OF 
REAUU0  /0|^NI 


OFF-HOOK  INPUT  CHANNELS  = ) 

740 

OF  OUTPUT  CHANNELS  =') 


IF(NI .EO.0)GO  TO 
WRITE (25,901)NI 
FORMAT(  + ,33X, 13/ '0NO. 

READ(10 .0ji  NO 
WRITE125/902) NO 

FoRMAT(  + ,25X,I3/  0ENTER  MESSAGE  SEGMENT  LIST:  ) 

READ(10,0) J 
WRITE (25,903) J 
FORMAT ( ,14) 

MSLST (I)=J 
1 = 1 + 1 

IFTj  .GT.0)iGO  TO  10 
WRITE(25,904) 

FORMAT  ( (30UTPUT  FILE  DESIRED  ?') 

CALL  YESNO(QFFLG) 

WRITEl25,9I3)  / 

FORMAT ( 0SEED  * ) 

READ(10,0) J 
IF(J.GT.0) ISEED»J 

WKITE{25,914) ISEED. 

FORMATC  + ,7X.I10/  01NITIALIZE 
CALL  YE§N^|^HlF) 

6sTA^GERED  PRIORITY  ROTATION 


CHANNEL  MODELS  ?') 


write (2 

FORMAT ( 

CALL  YESNO (ROTFG) 

wRlTEUS/Sl^l 

FORMAT ( 0ADAPT1VE 
CALL  YESNO (ADTFG) 


? ) 


THRESHOLD  ?') 
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1 


2 

905 

3 

4 

5 

6 

906 

7 

10 

11 

907 

12 

il 

908 

15 

16 
17 

917 

20 

21 

22 

911 

23 

24 

918 

25 

26 

27 

912 

30 

31 

13 

32 

990 

33 

34 

35 

909 

36 

37 

920 

40 

41 

42 

921 

43 

44 

780 

45 

46 

922 

50 

51 

52 

923 

53 

54 

55 

56 

57 

924 

60 

61 

62 

63 

925 

OF  UPDATE  REQUEST  DURING  SPURT  = ) 


wkITfc{25  905) 

FORMAT  (iScHANGE  OTHER  PARAMETERS  ?') 

CALL  YESNO(PFLG) 
if7.n.pflg)GO  to  13 
WRITE(2S,906) 

FORMAT ( 0MEAN  SPURT  LENGTH  IN  SECS  = ) 

READ(10 ,0) SPLNG 
WRITE (25,507) SPLNG 

FORMAT!  + ,28X,F5.2/  0SPEECH  ACTIVITY  = ) 

READ{10,0) ALPHA 
WRITE '25, 508IALPHA 
forma! ( + ,18X,F6.3) 

IF(ADTFG)GO  TO  13 
WRITE (25,917) 

FORMAT!  0THRESHOLD  = ) 

READ(1(3,0)  DELTA 
WR1TE(25,911)DELTA 
FORMAT!  + .12X,I3) 

WRITE (25,918) 

FORMAT!  ^PROBABILITY 
READ(10,0) PURDS 
WRITE (25,912)  PURDS 
FORMAT!  + ,45X,F5.2) 

WRITE!25,990) 

FORMAT!  0SPEECH  SIGNAL  GAIN  IN  DB  = ) 
READI10,0TspG 
WRITE (25,909) ^PQ 

FORMAT!  + ,27X,F6.1) 

WRITE (25,920) 

FORMAT!  0SIMULATE  SAW  BIT  ERRORS  ? ) 

CALL  YESNOlSBEFG) 

WA1TE(25  92l) 

FORMAT!  (Ss I MULATE  PCM  BIT  ERRORS  ?') 

CALL  YESNO(PBEFG) 

IFT.N.TsbEFG.O.PBEFG) ) GO  TO  20 
IF(CRDFG)  go  to  20 
WRITE (25,922) 

FORMAT!  0PROBABILITY  OF  ERROR  = ) 

READ  (10,0)1  PE 
WR1Tl(25  923)P£ 

FORMAT!  +,24X, Ell. 3/'0TEST  CHANNEL  POSITION  IN  SAW  =') 

READ(10 ,0) NTCPS 

IF(NTCPS.Lt.1)NTCPS=1 

IF(NTCPS.GT.48)NTCPS=48 

wRITE(25,924)NTCPS 

FORMAT!  + ,32X,I3/  0ESEED  = ) 

READ (10,0) 1 
IF(I.GT.0JESEED=I 
WRITE(25,525)ESEED 
FORMAT!  + ,8X,I10) 
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C INITIALIZE: 

C NT  = PRIORITY  OF  TEST  CHANNEL  (ROTATES) 

20  NT=1 

C NP  = CHANNEL  WITH  HIGHEST  PRIORITY  IN  CURRENT  FRAME 
NP  = NI 

C SET  UP  ADAPTIVE  THRESHOLD  RULE: 

IF(NO,LQ.16)GO  TO  389 
IfInO.E0.12  GO  TO  382 
C ASSUME  NO  = 18: 

DO  387  1=1,48 
ADTHR(I)=ADTHR18iI) 

387  KZ0A(I)=URNR*PURDA18 (I) 

GO  TO  381 
C NO  = 16: 

389  DO  385  1=1,48 

ADTHRID  =a6tHR16  (I) 

385  KZOA(I)=URNR*PURDAi6 (I) 

GO  TO  381 
C NO  = 12: 

382  DO  383  1=1,48 
AUTHR(I)=ADTriR12 (I) 

383  KZ0A(I)=URNR*PURDA12(I) 

381  CONTINUE 

C INCREMENT  FOR  STAGGERED  PRIORITY  ROTATION: 

DNT=NI-1 

DSGN=1 

CALL  OSDR(RBUF,MSLST) 

IF(OFFLG) CALL  OSDW(WBUF) 

C SAMPLE  COUNTER: 

N = 0 

C MEAN  SQUARE  SIGNAL  AND  ERROR  ACCUMULATORS: 

SQE=0 . 

SQS=0. 
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C NO.  OF  SAWS  IN  ERROR: 

NSERR=M 

C NO.  OF  DETECTED  SAW  ERRORS: 

NDSER=»3 

C NO.  OF  BAD  UPDATES: 

6 NBU=0 

7 C NO.  OF  TEST  CHANNEL  FREEZEOUTS  DUE  TO  SAW  ERRORS: 

It)  NFSE=0 

11  C NO.  OF  TEST  CHANNEL  SAMPLES  WITH  ERRORS: 

12  NPERR=0 

13  C INITIALIZE  BIT  ERROR  GENERATOR: 

14  AURNR=URNR/.9 

15  BURNR=AL0G( .9)/PE 

16  NERL=0 

17  CALL  EGEN 

20  C NO.  OF  BITS  PROCESSED  BY  BIT  ERROR  GENERATOR: 

21  NBCNT=0 

22  C NO.  OF  BIT  ERRORS: 

23  NECNT=0 
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C COUNTERS  OF  TEST  CHANNEL  FREEZEOUTS  AND  POSSIBLE  FREEZEOUTS: 
NPFO=0 
NFO=0 

C ACTIVITY  (NO.  OF  VOICE  SWITCHES  ON)  HISTOGRAM: 

DO  22  1=1,48 
22  ACTIV(I)=0 

IF(PFLG  CALL  QSCAL(QSCL) 

c MEAN  Spurt  length  in  Samples: 

SSPL=SPLNG/DT 

C MEAN  GAP  LENGTH  IN  SAMPLES: 

^ SGPL= (1. -ALPHA) *SSPL/ALPHA 

C THRESHOLD  TO  GENERATE  ZERO-ONE  RANDOM  VARIABLE 
C WITH  P(1)=ALPHA: 

KZO=URNR*ALPHA 

C NUMBER  OF  SIMULATED  CHANNELS: 

NSC=NI-1 

C START  SIMULATED  CHANNELS  WITH  STEADY  STATE  DISTRIBUTION: 

IF] .N.CHIF) GO  TO  27 

C ACTIVITY  COUNT  - START  AT  1 SINCE  TEST  CHANNEL  ALWAYS  ACTIVE* 
NACTV=1 
DO  25  1=1, NSC 
CALL  RAND(ISEED,J) 

IF] J .GT.KZO)GO  TO  26 
SrATL(I)=.T. 

CNT(I) =EXPRN ( .T.) 

NACTV=NACTV+1 
GO  TO  25 

26  STATE(I)=.F. 

CNT(1)=EXPRN(.F.) 

25  CONTINUE 

C TEST  CHANNEL  SAMPLE  IN  ENCODER  & DECODER  FRAME  MEMORIES: 

27  LTCUD»128 
LhCOD=128 

C CONSTANT  FOR  UPDATE  REQUEST  0-1  RV: 

KZO-URNR*PURDS 
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C MAIN  LOOP: 

C GET  SPEECH  SAMPLE: 

50  CALL  SGET(J.END) 

C IF  NONE  LEFT,  EXIT  LOOP: 

IF(END)GO  TO  200 
C COUNT  SAMPLES: 

N = N + 1 

SAMP=J 

FSAMP=SAMP 


C ATTENUATE  SIG  BY  -SPG  DB 

FSAMP=FSAMP*10 (SPG/20 . ) 

C FMINl  ADDED  SO  THE  A CAIN  ] 0 WILL  NOT  CAUSE  FRACTIONAL 
C OVERFLOW  WHEN  FSAMP  IS  CONVERTED  TO  SAMP 
FSAMP=AMAXl (-1 . .FSAMP) 

SAMP=FMIN1 (1 . .FSAMP) 

C ACCUMULATE  SICSNAL  POWER  IN  TEST  CHANNEL: 

SQS=SQS+FSAMP*FSAMP 
C QUANTIZE  TEST  CHANNEL  SAMPLE: 

PCODE=PCMWD (SAMP) 

C FORM  MAGNITUDE  DIFFERENCE  BETWEEN  PRESENT  AND  SAVED  CODES: 
DCODE=IABS (PCODE-LTCOD) 

C IF  THRESHOLD  IS  ADAPTIVE,  SET  THRESHOLD  BASED  UPON  CURRENT 
C ACTIVITY  & UPDATE  REQUEST  CONSTANT  BASED  UPON  THRESHOLD: 
IF( .N.ADTFG)GO  TO  53 
DELTA=ADTHR (NACTV) 


KZO=KZOA (NACTV) 

IF  CURRENT  PRIORITY  OF  TEST  CHANNEL  IS  1 THRU  NO.  OF  OUTPUT 
CHANNELS,  IT  WILL  BE  UPDATED  IF  IT  WANTS  TO  BE: 

IF (NT. LE. NO) GO  TO  70 

OTHERWISE,  SET  UP  TO  COUNT  SIMULATED  CHANNELS  NEEDING  UPDATE: 
I=NP 


C UPDATE  REQUEST  COUNTER: 
J=0 


c VOICE  Switch  of  channel  i on  ? 

55  IF( .N.STATE(I) )GO  TO  57 

C IF  SO.  GENERATE  0-1  RV  FOR  UPDATE  REQUEST: 
CALL  RAND(ISEED,K) 

1F(K. LT.KZO) J=J+1 
C ALL  OUTPUT  CHANNELS  USED  YET  ? 

IF(J.GE.NO)GO  TO  60 
C NO.  GO  TO  NEXT  CHANNEL: 

57  1=1+1 


C IS  THIS  NEXT  CHANNEL  THE  TEST  CHANNEL  ? 
IF(I.LT.NI)GO  TO  55 
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70 
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7 

10 
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11 

12 

c 

13 

14 

c 

15 

16 

17 

c 

20 

21 

c 

22 

23 

c 

24 

c 

25 

26 

27 

205 

30 

31 

C 

32 

210 

V 

C 

34 

35 

36 

37 

40 

C 

41 

260 

42 

YES.  THERE  IS  AN  OUTPUT  SLOT  FOR  THE  TEST  CHANNEL 
IF(DCODE. LE.DELTA)GO  TO  65 
UPDATE  TEST  CHANNEL  IN  ENCODER  FRAME  MEMORY  IF 
THRESHOLD  wAS  EXCEEDED: 

LTCOD=PCODE 

SAW  BIT  ERRORS  BEING  SIMULATED  ? 

IF( .N.SBEFG)GO  TO  230 
IF  SO,  GENERATE  SAW  ERROR  PATTERN 
CALL  SAwEP 
ANY  ERRORS  ? 

IF( .N.SERRF)GO  TO  230 
YES.  SIMULATE  SAW,  ETC. 

TCSB=1 
CALL  SAWGN 

DETECTABLE  ERROR  ? (CORRECT  NUMBER  OF  1 S) 
IF(NRSl.NE.NO)GO  TO  260 
NO.  TEST  CHANNEL  BIT  IN  ERROR  ? 

IF(ESAW(NTCPS) .EQ.0)GO  TO  260 
NO.  POSITION  OF  TEST  CHANNEL  BIT  IN  SEQUENCE  OF  1 
THE  SAME  IN  TRANSMITTED  & RECEIVED  SAWS  ? 

DO  205  1I=1,NTCPS 
KK=KK+SAW ( Ill-ESAW ( II) 

IF(KK.EQ.0)GO  TO  230 
NO.  BAD  UPDATE: 

NBU=NBU+1 

GENERATE  RANDOM  SAMPLE: 

CALL  RAND(ESEED,KK) 

L.  COD=KK  . A.  377B 
GO  TO  280 

r-RLEZEOUT  DUE  TO  SAW  ERROR: 

NFSE=HFSE+1 
GO  TO  280 


S 


A-9 


SFLC3 


VEH  1 


RE\/  A 1 OCT  76 


PAGE  11 


1 

C 

2 

230 
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11 

12 

13 

14 

233 

15 

C 

16 

17 

C 

20 

21 

C 

22 

23 

24 

25 

C 

26 

240 

27 

30 

31 

C 

32 

C 

33 

60 

34 

C 

35 

36 

C 

37 

65 

40 

C 

41 

42 

C 

43 

44 

C 

45 

46 

47 

C 

50 

51 

C 

52 

PCM  BIT  ERRORS  BEING  SIMULATED  ? 

IF( .N,PBEFG)GO  TO  240 
ytS.  GENERATE  8 BIT  ERROR  PATTERN: 

PCMEP=0 
PERRF=. F. 

DO  233  11=1,8 
PCMEP= . B. PCMEP 
CALL  EGEN 

1F(EBIT.EQ.0)GO  TO  233 
PERRF=.T. 

PCMEP=PCMEP.O. 1 
CONTINUE 
ANY  ERRORS  ? 

1F( .N.PERRF)GO  TO  240 
YES.  COUNT  PCM  ERROR  ON  TEST  CHANNEL: 

NPERR=NPERR+1 

PUT  CORRUPTED  SAMPLE  IN  DECODER  FRAME  MEMORY: 

LRCOD=LTCOD . XOR . PCMEP 
GO  TO  280 

ERROR  FREE  PCM  SAMPLE  GOES  TO  DECODER: 

LRCOD=LTCOD 
GO  TO  280 

ALL  OUTPUT  CHANNELS  FILLED  BEFORE  TEST  CHANNEL  WAS  REACHED. 
COUNT  POSSIBLE  FREEZEOUT: 

NPFO=NPFO+l 

IF  TEST  CHANNEL  WANTS  UPDATE,  COUNT  A FREEZEOUT: 

IF(DCODE.GT. DELTA) NFO=NFO+1 
SAW  BIT  ERRORS  BEING  SIMULATED  ? 

1F{.N.SBEFG)G0  TO  280 
YES.  GENERATE  SAW  ERROR  PATTERN: 

CALL  SAWEP 

any  errors  ? 

IF( .N.SERRF) GO  TO  280 
YES.  SIMULATE  A SAW,  ETC. 

TCSB=0 
CALL  SAWGN 
DETECTABLE  ERROR  ? 

1F(NRS1.NE.N0)G0  TO  280 

NO.  IF  TEST  CHANNEL  BIT  IS  IN  ERROR,  DO  A BAD  UPDATE: 
IF(ESAW(NTCPS) .NE.0)GO  TO  210 
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C RECONSTRUCT  TEST  CHANNEL  AT  SPEC  DECODER: 

280  QbMP=PCMDC (LRCOD) 

C ACCUMULATE  ERROR  POWER 
X=FSAMP-QSMP 

sul=sqe+x*x 

C OUTPUT  SAMPLE  TO  DISK  IF  DESIRED: 

IF{OFFLG) CALL  SPUT(QSMP) 

C ROTATE  PRIORITY: 

IF(ROTFG)GO  TO  81 
C STANDARD  ROTATION  SCHEME: 

Nr=NT+l 

IF(NT.GT.NI)NT=1 
GO  TO  82 

C STAGGERED  ROTATION  SCHEME: 

81  NT=NT+DSGN*DNT 
DNT=DNT-1 
DSGN=-DSGN 

IF(DNT.GE.0)GO  TO  82 
NT=1 

DNT=NI-1 

DSGN=1 

82  NP=NI+1-NT 

C ACCUMULATE  ACTIVITY  HISTOGRAM: 

ACTIV (NACTV) =ACTIV (NACTV) +1 
C UPDATE  CHANNEL  MODELS: 

DO  90  1 = 1, NSC 

C DECREMENT  SPURT/GAP  LENGTH  COUNTER: 

CNT(1)=CNT(I)-1 
C END  OF  SPURT/GAP  ? 

IF(CNT(I) .GT.0)GO  TO  90 
C YES.  SWITCH  TO  OPPOSITE  STATE: 

STATE ( I ) = . N , STATE ( I ) 

C GENERATE  EXPONENTIALLY  DISTRIBUTED  RV  FOR  LENGTH  OF 
C NEXT  SPURT/GAP: 

CNT(I)=EXPRN(STATE(I) ) 

C UPDATE  NO.  OF  CHANNELS  ACTIVE: 

NACTV=NACTV-1 
IF (STATE (I) ) NACTV=NACTV+2 
90  CONTINUE 
C LOOP: 

GO  TO  50 
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C END  OF  TEST  CHANNEL  INPUT  FILE: 

200  IF (OFFLG) CALL  CSDW 

SNR=10 . *ALOG10 (SQS/SQE) 

SQS=SOS/N 

C LBSIG  AND  DbSIG2  ADDED  9/14/76  TRM 

C DB5IG2  IS  THE  DB  OF  THE  SIG  RELATIVE  TO  A FULL  WAVE  SINUSOID 
DBSIG=10. *ALOG10 (SQS) 

UBSIG2=DBSIG+3 . 

X=100 . *NFO/N 

Y=100 . ‘NPFO/N 

WHITE (25,950) N, SNR, NFO,X,NPFO,Y 
WHITE  25,952) ACTIV 
WHITE (25,955) ISEED 

WHITE (25.962) DBSIG,DB§IG2 

962  FORMAT ( 0SIGNAL  PWR  = ,F6.1.  DB  ,/,  0§IGNAL  PWR  RELATIVE 
1,  TO  A FULL  WAVE  SINUSOID  = ,F6.1,  DB  ) 

ISEED=7654321 
IF(CRUFG)GO  TO  770 
GO  TO  760 

740  IF ( .N.  (SBEFG.O.PBEFG) ) CALL  EXIT 
X=NBCNT 
X=NECNT/X 

WHITE  (25,960) NBCNT,NECNT,X ,NPERR,NSERR,NDSER,NBU,NFSE , ESEED 
CALL  EXIT 

960  FORMAT('0NO.  OF  ERROR  PATTERN  BITS  GENERATED  =',110/ 

1 , NO.  OF  BIT  ERRORS  = ,110/ 

2 , ACTUAL  PROBABILITY  OF  ERROR  = ,E12.3/ 

3 , NO.  OF  TEST  CHANNEL  SAMPLES  WITH  BIT  ERRORS  = ,110/ 

4 , NO.  OF  SAWS  IN  ERROR  = *110/. 

5 , NO.  OF  DETECTED  SAW  ERRORS  = .110/ 

6 , NO.  OF  BAD  TEST  CHANNEL  UPDATES  = ,110/ 

7 , NO.  OF  TEST  CHANNEL  FREEZEQUTS  DUE  TO  SAW  ERRORS  = ,110/ 

8 LAST  SEED  FOR  BIT  ERRORS  = ,110) 

950  FgRMAT( '0  ' ,110, ' SAMPLES '/' 0S/N  =',F7.2,'  DB'/ 

1 ,0NUMBER  OF  TEST  CHANNEL  FREEZEOUTS  = ,110/ 

2 , PERCENT  TEST  CHANNEL  FREEZEOUTS  = ,F7.2/ 

3 ,0NUMBER  OF  POSSIBLE  TEST  CHANNEL  FREEZEOUTS  = ,110/ 

4 PERC^iNT  POSSIBLE  FREEZEOUTS  = ,F7.2) 

952  FORMATT  0ACTIVITY  HISTOGRAM;  /(  ,6110)) 

955  format!  0LAST  SEED  = ,111) 
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C EXPONENTIAL  RANDOM  NUMBER  GENERATOR 
INTEGER  FUNCTION  EXPKN(SPFG) 
LOGICAL  SPFG 
CALL  RAND(ISEED,JJ) 

C TRUNCATE  TO  10* MEAN; 

IF(JJ .LT.6093) JJ=6093 
XX=-ALOG (JJ/URNR) 

C SELECT  ONE  OF  TWO  MEANS: 

IF (SPFG) GO  TO  850 

EXPRN=SGPL*XX 

RETURN 

850  EXPRN=SSPL*XX 
RETURN 


SUBROUTINE  YESNO(FLAG) 
LOGICAL  FLAG 
800  CONTINUE 

A JPSR,  SICC 

A ARMD  L:  CHAR:  0 

IF(CHAR.EQ. 1RY)G0  TO  810 
IF (CHAR. EQ. IRN) GO  TO  820 
WRITE (25,899) 

GO  TO  8i40  , 

899  FORMAT ( ? ) 

810  FLAG=.T. 

WRITE (25,898) 

RETURN  , 

898  FORMAT ( YES  ) 

820  FLAG=.F. 

WRITE (25,897) 

return  , 

897  FORMAT ( NO  ) 
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C BIT  ERROR  GENERATOR: 

SUBROUTINE  EGEN 
C COUNT  BITS  PROCESSED: 

NBCNT=NBCNT+1 

C DECkEMENT  run  length  to  next  ERROR: 

NERL=NERL-1 

EBIT=k) 

IF(NERL.GE.0) RETURN 
C TIME  FOR  AN  ERROR: 

EBIT=1 

NECNT=NECNT+1 

C GENERATE  EXPONENTIALLY  DISTRIBUTED  RUN  LENGTH 
C TO  NEXT  ERROR: 

CALL  RAND(ESEED, JJJ) 

C TRUNCATE  TO  PREVENT  VERY  LONG  RUNS  & INTEGER  OVERFLOW: 
IFCJJJ .LT.5484) JJJ=5484 

C NOTE  SCALING  TO  AVOID  ARGUMENTS  OF  ALOG  NEAR  1. 

C AND  POSSIBLE  ROUNDOFF  ERROR  ON  SHORT  RUNS  TO 
C NEXT  ERROR: 

NERL=BURNR-ALOG ( JJ J/AURNR) /PE 
RETURN 


C 

c 

c 

c 


4k)0 

C 


SArt  ERROR  PATTERN  GENERATOR: 

SUBROUTINE  SAWEP 
CLEAR  SAW  ERROR  FLAG: 

SERRF=.F. 

GENERATE  48  BIT  ERROR  PATTERN  & SET  FLAG 
IF  ANY  ERRORS: 

DO  400  111=1,48 
CALL  EGEN 

IF (EB1T.NE.0) SERRF=.T. 

ESAW (III) =EB1T 
COUNT  SAW  ERRORS: 

IF(SERRF) NSERR=NSERR+1 
RETURN 
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C SAW  SIMULATION: 

SUBROUTINE  SAWGN 
C CLEAR  OLD  SAW: 

DO  410  111=1,48 
410  SAW(I1I)=0 

C PUT  A 1 IN  TEST  CHANNEL  POSITION  SO  NO  RANDOM 
C BITS  WILL  BE  PLACED  THERE: 

SAW (NTCPS) =1 

C NUMBER  OF  RANDOM  BITS: 

JJ=NO-TCSB 

C GENERATE  THIS  NUMBER  OF  BITS  iN  RANDOM  POSITIONS: 

DO  415  III=1,JJ 

C GENERATE  A RANDOM  NUBER  BETWEEN  1 5.  48: 

416  CALL  RAND(ESEEU,KK) 

KK= ( (KK.R.3) .A.77B  +1 
IF(KK.GT.48jGO  TO  416 

C IF  THERE  IS  ALREADY  A BIT  THERE,  DO  IT  AGAIN 
IF(SAW(KK) .NE.0)GO  TO  416 
415  SAw(KK)=l 

C SLT  THE  TEST  CHANNEL  BIT  TO  WHAT  IT  SHOULD  BE: 

SAW ( NTC  PS ) =TCSB 

C FORM  THE  RECEIVED  SAW  & COUNT  THE  NUMBER  OF  1 'S  IN  IT: 
NRS1=0 

DO  420  11=1,48 
KK=SAW(II) .XOR. ESAW (II) 

NKS1=NRS1+KK 
420  ESAW(II)=KK 

C COUNT  A DETECTED  ERROR  IF  NUMBER  OF  1 'S  IS  WRONG: 
IF(NRSl.NE.NO)NDSER=NDSER+l 
RETURN 

end 
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'I 

J 

4 

5 

6 
7 

Ik) 

11 

12 

li 

14 

15 
lb 
17 

2k) 

21 

22 

24 

25 

26 
27 
Jk) 

32 

33 

34 
3 5 

36 

37 
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41 

42 

43 

44 

45 

46 

47 
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62 
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64 
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66 

SI 

71 

72 

73 


P^OGRArt  MCDAQ2 

C RiULTICHANNEL  DIFFERENTIAL  ADAPTIVE  QUANTIZATION 
vARIAtiLL  WORD  LENGTH  ADPCM 

c uasld  on  jayant  s algorithm 

C WORD  length  is  varied  BASED  ON  ACTIVITY,  I.E., 

C NUMBER  OF  VOICE  SWITCHES  ON. 

C ACTIVITY  IS  SIMULATED  AS  IN  SPECl  PROGRAM. 


C 2 TO  8 BITS  PER  SAMPLE 
C FIXED  RATIO  STEP  SIZES: 

L 64  STEP  SIZES  wlTH  A RATIO  OF  2**. 125 

C ADAPTATION  RULE  IS  APPROXIMATELY  THAT  GIVEN  BY  JAYANT 
C FOR  A 5 BIT  SYSTEM  WHEN  ThE  WORD  LENGTH  IS  5 OR  GREATER. 

C FOh  smaller  WORD  LENGTHS,  THIS  STEP  SIZE  ADJUSTMENT  TABLE 
C IS  SUBSAMPLED. 


C R.W.S.  7/22/76 
C VERS.  1,  kEV.  a 


C 


FRACTION  SAMP , FR , FSD , FRL , FD , DF 
DATA  DF/IB/ 

INTEGER  RBUF(1664:0) ,WBUF(1664) ,MSLST(20) 
LOGICAL  END 

INTEGER  IRTBL{64:0) ,IDCNT(16:0) 

COUNTER  INCREMENTS: 

DATA  IDCNT  /8 *-l , 2 , 5 , 7 , 9 , 1 0 , 1 1 , 1 3 , 1 4/ 


LOGICAL  STATEa00)  , CH  IF  , NWST , CRDRD  , MORE  , OFI LE 
INTEGER  CNT(100)  ,WLH(8)  ,EXPRN,A(20)  ,B(20)  ,C(20) 
DIMENSION  D(20)  ,E  (20)  ,Fj20) 

DATA  URNR/134217728 ./ ,DT/125. E-6/ 

ISEED=7654321 


WRITE (25f 800) 

800  FORMAT ( MULTIPLE  INPUT?  ) 
CALL  YESNO (CRDRD) 

IF(.N. CRDRD)  GO  TO  810 


Bdl 
88  5 

883 


OFILE=.F. 
1 = 1 


WRITE (25,881) 

FuRMAT(  SEG, INPUT, BITS, SPURT, ACTIV, GAIN 
READ(10,k)T  A(I),B(l),C(ir,D(l[,E(I),F(I) 
WRITE (25,&83)  A(I),B(I),C(I),D(I),E(I),P 


1=1+1  , , 

FORMAT ( 0 ,3I5,3F10.5, 
CALL  YESNO (MORE) 

IF (MORE)  GO  TO  885 
MAX=I-1 


MORE  ?',//) 


??? 

(I) 


,//) 


860 


12*1 

1F(I2.GT.MAX) 

MSLST(1)»A(I2) 

Nl*B(ll 

NBPF»C(I2) 

SPLNG»D(I2) 

ALPHA«Ei 12) 


GO  TO  850 


A-16 


MCDAQZ  VER 


REV  A 18  OCT  76 


7 4 

75 

76 
7 7 

1 uu 
ikJl 

ik)3 

lk)4 

105 


110 

111 

112 

113 

114 

115 
11b 
117 
120 
121 
122 

123 

124 

125 

126 
127 


SPG=f (12) 

ISEED=7654321 

12=12+1 

rtklTE  (2^,830)  MSLSTd)  , N1 , NBPF  . SPLNG,  ALPHA , SPG 
830  FORMAT ( ISEG  NI  SPEECH  BITS  BURST  ACTIV 
1/,1H0,I3,2X,13,6X,I3,8X,3F7.3) 

CHIF= .T. 

GO  TO  840 


810 

792 

903 
25 

904 


WRITE (25/792) 

FoRMAT(  0OUrPUT  FILE??') 

CALL  VESNu(OFILE) 

WRITE (25,903) 

FORMAT ( 0ENTER  MESSAGE  SEGMENT  LIST:  ) 

READ (10,0) J 
WRITE (25,904)3 
FORMAT ( ',14) 

MSLST(I)=J 
1 = 1 + 1 

IF(J.GT.0)GO  TO  25 


c Initialize  sample  data  i/o: 

840  CALL  OSDR(kBUF,MSLST) 
IF(OFILE)  GO  TO  842 
CALL  OSDW(WBUF) 

842  IF(CRDRD)  GO  TO  888 


> 
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1 

2 

i 

4 

6 

6 

7 

m 
11 
12 
I J 

14 

15 

16 
17 
26 
21 
22 

24 

2d 

26 

27 

jk) 

31 

32 

33 

34 

35 

36 

37 
46 

41 

42 

43 

44 

45 

46 

47 


VvRlTE  (25,920) 

920  FORMAT ( 0NO.  OF  OFF-HOOK  INPUT  CHANNELS  =') 
R£AD(10.0)M 

WRITE  (25,^211,NI  ^ . 

921  FORMAT!  + ,33X,I3/  0IN1TIALIZE  CHANNEL  MODELS  ? ) 
CALL  YESNO(CHIF) 

kvRITE]25,922) 

922  FORMAT!  0SEED  = ) 

READ(10,W) I 

IF(I .GT.0) ISEED=1 
wR1TE(25.323) ISEED, 

923  FORMAT!  + ,7X,110/  kJMEAN  SPURT  LENGTH  IN  SECS.  =' 
READ!10.e) SPLNG 

wRlTE(25,924)SPLNG 

924  FORMAT!  + ,28X,F5.2/  0SPEECH  ACTIVITY  =') 
RLAU(1O.0) ALPHA 

vikITE  (25.925)  ALPHA 

925  FORMAT!  + ,18X,F6.3/ 

1 0WO.  OF  HITS  PER  FRAME  = ) 

REAb(10,0lNbPF 
wRITE!2d,926)NBPF 

926  FORMAT!  + ,25X,15,// 

READ(10 .0) SPG 
wRITLl25.y27)  SPG 

927  FORMAT!  0GAIN=  ,F7.2,  DB  ,//) 


'0SPEECH  SIGNAL  GAIN  IN  DB? 


866  CONTINUE 

C INITIALIZE  irvORD  LENGTH  HISTOGRAM: 

DO  22  1=1,8 
22  WLH!I)=0 


C MEAN  SPURT  LENGTH  IN  SAMPLES: 

3SPL=SPLNG/DT 

C MEAN  GAP  LENGTH  IN  SAMPLES: 

SGPL=(1 .-ALPHA) *SSPL/ALPHA 
C NO.  OF  SIMULATED  CHANNELS: 

NSC=N1-1 

C THRESHOLD  TO  GENERATE  0-1  RV  WITH  P!1)=ALPHA: 
KZO=URNR*ALPHA 
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1 

C 

z 

3 

c 

4 

5 

6 

7 

10 

11 

12 

13 

14 

31 

16 

30 

17 

20 

C 

21 

40 

22 

23 

25 

36 

26 

27 

C 

51 

c 

32 

33 

54 

c 

35 

36 

37 

40 

41 

c 

42 

c 

43 

50 

44 

45 

C 

46 

47 

50 

51 

52 

53 

54 

C 

55 

START  SIMULATED  CHANNELS  wITH  STEADY  STATE  DISTRIBUTION: 

IF( .N.CHIF)GO  TO  40 

ACTIVITY  COUNT  - START  AT  1 SINCE  TEST  CHANNEL  ALWAYS  ACTIVE 
NACTV=1 


DO  30  1=1, NSC 
CALL  RAND(ISEED,J) 
IF (J .GT.KZO)GO  TO 
STATE (I) = .T. 

CNT (I) =EXPRN ( .T. ) 
NACTV=NACTV+1 

GO  TO  3 k) 

STATE ( I ) = . F . 

CNT (I ) =EXPRN ( .F. ) 
CONTINUE 


31 


BUILD  QUANTIZER  RANGE  LOOKUP  TABLE: 
DX=2.** .125 
X=16383. 

DO  36  l = 63,k),-l 
IRTBL ( I ) =X 
X=X/DX 


INITIALIZE  SAMPLE  COUNTER: 

N = k) 

INITIALIZE  MEAN  SQUARE  SIGNAL  S.  ERROR  ACCUMULATORS: 
SQS  = ij  . 

SQE=6 . 

Initialize  prediction  & quantizer  range: 

IY  = k) 

ICNT=U 

FR=IRTBL(ICNT) 

MAIN  LOOP: 

GET  NEXT  SAMPLE,  IF  THERE  IS  ONE: 

CALL  SGET{IX,END) 

IF(END)GO  TO  200 
ATTENUATE  SIGNAL  BY  -SPG 
XXX=IX 

XXX=XXX*10. ** (SPG/20.) 

IX=XXX 

IIX=IX 

IFdX.GT. 16384)  1X=16384 
IFdX.LT. -16334)  IX=-16384 
COUNT  SAMPLES: 

N = N + 1 
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1 

2 

3 

4 

D 

6 

7 

1 

1 

12 

13 

14 

15 

16 
17 
20 
21 
22 

23 

24 
2o 
26 
27 

30 

31 

32 

33 

34 
3d 
36 
3 7 
4ti 


C SET  wOkD  LENGTH  BASED  ON  ACTIVITY: 

IBPS=NBPF/NACTV 

C Don  T allow  MORE  THAN  8 BITS/SAMPLE: 

IF(IBPS.GT.8) 1BPS=8 
C ACCUMULATE  WORD  LENGTH  HISTOGRAM: 

WLH (IBPS)=WLH (IBPS)+1 
C SET  UP  ADPCM  PARAMETERS: 

ML=1.L.IDPS 

IM1D=ML.R. 1 

ML=ML-1 

IRSHF=15-IBPS 

ILSriF=IRSHF-l 

lBPS5=IbPS-5 

NBPS5=-IBPS5 

C FORM  PREDICTION  ERROR: 

ID=IX-1Y 
C gUANTIZE  IT: 

FD=1D 

C CLAMP  TO  QUANTIZER  RANGE: 

IF (FU.GT. FR) FD=FR-DF 

1F(FD.LT.-FR)FD=-FR+DF 

FSD=FD/FR 

IQ= ( ( (FSD.R. IRSHF) .A.ML)+IMID) .A. ML 
C FORM  representation  LEVEL: 

FHL= ( ( .B. (IQ-IHID) ) +1) .L. ILSHF 
FRL=FRL*FR 

C RECONSTRUCTED  SAMPLE  = REPRESENTATION  LEVEL  + PREDICTION: 
1XH=FRL 
IXri=lXri+l Y 

o NEXT  PREDICTION  = .875  * RECONSTRUCTED  SAMPLE: 

IY=IXH- (IXH.R.3) 


C>Lr  LP  UjU^U;OjNJ  NlN-  tOh-'h-H-h-'h-  H'H-'h- '-OO^  CP•t^  U.  N.'»— 

S-oa  CniJfcU>K‘H-S  ^OO'LPi^U  hJH-K  MH-C  LP*^  U' f^*  S. 
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S.. 

c 


c 

63 

C 

C 

61 


UPDATE  STcP  SIZE  FOR  NEXT  TIME  BASED  ON  CURRENT  LEVEL: 
FORM  MAGNITUDE  OF  LEVEL  STARTING  AT  0: 


IQM=IQ 

iF{ir 


(IQM. LT. IMID) IQM=ML-IQM 
IQW=IQM-IMID 
IF(IBPS5) 61,62,63 

WORD  LENGTH  . GT . 5 BITS  - KEEP  4 MSB  S 
IQM=IOM. R. IBPS5 
GO  TO  62 

WORD  LENGTH  .LT.  5 BITS  - POSITION  IN  UPPER 
PART  OF  4 BIT  WORD: 

IQM=IOM. L.NBPS5 


62  ICNT=ICNT+IDCNT(IQH) 
IF(ICNT.LT.0) ICNT=0 
IF(ICNT.GT.63) ICNT=63 
FR=IRTBL(ICNT) 


C OUTPUT  SAMPLE  TO  DISK: 

IF(.N.OFILE)  go  to  8000 
CALL  SPUT(IXH) 

C ACCUMULATE  SIGNAL  & ERROR  POWERS: 
80U0  SAHP=IIX 
X=SAMP 
SCS=SQS+X*X 
SAMP=IIX-IXH 
X=3AMP 
SQE=SQE+X*X 

C UPDATE  CHANNEL  MODELS: 

DO  90  1=1, NSC 

C DECREMENT  SPURT/GAP  LENGTH  COUNTER: 

CNT(I)=CNT(I)-1 
C END  OF  SPURT/GAP  ? 

IF(CNT(I) .GT.0)GO  TO  90 

c iES.  Switch  to  opposite  state: 

NWST=.N.STATE (I) 

STATE ( I ) =NWST 

C GENERATE  LENGTH  OF  NEXT  SPURT/GAP: 

CNT(I)=EXPRN (NWST) 

C UPDATE  ACTIVE  CHANNEL  COUNT: 
NACTV=NACTV-1 
IF (NwST) NACTV=NACTV+2 
90  CONTINUE 


C LOOP: 

GO  TO  50 
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0 

*1 

12 

13 

14 

15 


200 


850 

950 


End  of  input  FILE: 

CALL  CSDw 
SNR=10.*ALOG10 (SC 
SPWR=10.*ALOG10 U 
WRITE (25,950)1 
WRITE (25,951)  : 

1F(CRDRU)G0  TO  860 

CALL  EXIT,  , , , 

FORMAT ( 0 ,110,  SAMPLES  / 0S/N  =' 
’ 0SIGN^L  POkyER  RELATIVE  TO  A FULL 


951 


^S/SOE) 

SQS/N)+3 

NR,SPWR, 


(1,WLH(I) 


2,Eb.2,  DB  ,///, 

3 0WORD, LENGTH  HISTOGRAM:  /( 
FORMAT ( 0LAST  SEED  = ,111//) 


,12 
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1=1,8) 

,F7.2,'  DB'/ 
WAVE  SINUSOID= 

,110)  ) 
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1 

/ 

8 

4 

5 

6 
7 

Ih] 

11 

12 

13 

14 

15 

16 
17 

2k) 

21 

22 

23 

24 

25 

26 
27 
3i) 
31 
3z 

33 

34 

35 

36 

37 
4ti 

41 

42 
41 


L EXPONENTIAL  RANUOKt  .<UrtBER  GENERATOR 

Integer  function  expkn(spfg) 

LOGICAL  3PFG 

CALL  RAND(ISEEL,JJ) 

C TRUNCATE  TO  1«)*NEAN: 

IF  (JJ  . LT.6k)93)  JJ=66  93 
XX=-ALOG(JJ/URNR) 

C SELECT  ONE  OF  TrtO  MEANS: 

1F(SPFG)G0  TO  85U 

EXPRN=SGPL*XX 

RETURN 

85i()  EXPRN  = SSPL*XX 
RETURN 


SUBROUTINE  YESNO(FLAG) 
LOGICAL  FLAG 
80k3  CONTINUE 

A JPSR,  SICC 

A ARMD  L;  CHAR:  0 

IF (CHAR. EG. IRY) GO  TO  810 
IF (CHAR. EQ. 1RN)G0  TO  820 
WRITE (25,899) 

GO  TO  800  , 

899  Format ( ? ) 

81k)  FLAG=.T. 

WRITE (25,898) 

RETURN  , 

898  FORMAT ( YES  ) 

820  FLAG= . F . 

write (25,897) 

RETURN  , 

897  tORMAT(  NO  ) 

END 


A-23 


APPENDIX  B 


# 


Seme  Details  of  the  SPEC  Slnulatlon 

1*  Stappered  Priority.  The  SPEC  3 program  simulated  systems 
which  agreed  In  every  respect  except  two  with  the  system  out- 
lined In  [1]». 

The  first  exception  Is  the  priority  arrangement.  In  [1] 
the  channel  priority  arrangement  assigned  the  highest  priority 
sequentially  to  channels  1,2,3# the  first  frame 
channel  1 has  highest  priority,  for  the  next  frame  channel  2 has 
the  highest  priority  and  so  on  (see  Pig.  5 of  [1]).  We  can  show 
that  the  overload  noise  power  can  be  reduced  If  we  use  the 
staggered  priority  rotation  scheme  1 , 6^ , 2 , £3, 3 , 62, • • • , There- 
fore this  priority  scheme  was  used  in  SPEC  3. 

2.  Threshold  settings.  The  second  way  In  which  the  SPEC  3 
simulations  varied  from  [1]  Is  In  the  SPEC  predictor  thresholds. 
The  threshold  Is  defined  by  a function  a(m)  where  a Is  the 
threshold  and  m Is  the  number  of  voice  switches  on  (l.e. , m is  the 
number  of  channels  which  request  an  update).  If  o Is  0 then 
every  channel  requesting  an  update  gets  one.  If  o'l  only  those 
channels  whose  current  sample  differs  from  Its  predicted  value 
by  more  than  1 gets  an  update.  In  general,  only  those  channels 
whose  current  sample  differs  from  Its  predicted  value  by  more 
than  a are  assigned  a transmission  slot. 

The  following  table  defines  the  function  a(m)  from  the 

SPEC  3 program.  Also  shown  In  the  table  are  the  probabilities  of 

exceeding  the  threshold  given  that  the  voice  switch  Is  on. 

•understanding  the  next  thj-ee  paragraphs  requires  a deeper  under- 
standing of  the  SPEC  system  than  Is  furnished  In  this  report. 
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The  quantities  In  this  table  were  used  In  simulating,  traffic 
In  the  SPEC  system. 


E 


Table  Bl.  Statistics  used  In  generating  traffic 
In  the  SPEC  3 program. 


No.  of  voice 
switches  on  n 

‘tS/o 

Threshold 

ue/2 

a 

i*2/6 

Prob.  of  e:tceed 
threshold  given 
voice  switch  Is 

148/0  146/2 

ine 

that 

on 

142/6 

1-12 

0 

0 

0 

.87 

.87 

.87 

13 

0 

0 

1 

.87 

.87 

.76 

lit 

0 

0 

2 

.87 

.87 

.68 

15 

0 

0 

3 

.87 

.87 

.69 

15 

0 

0 

3 

.87 

.87 

.59 

17 

0 

1 

u 

.87 

.76 

.53 

18 

0 

1 

14 

.87 

.76 

.53 

19 

1 

2 

14 

.76 

.66 

.53 

20 

3 

14 

.76 

.59 

.148 

21 

2 

3 

5 

.68 

.58 

.148 

22 

2 

3 

5 

.68 

.58 

.14 14 

23 

3 

3 

6 

.59 

.58 

.1414 

2U 

3 

k 

6 

.59 

.53 

.1414 

25 

7 

k 

6 

.59 

.53 

.1414 

26 

3 

k 

6 

.59 

.53 

27 

U 

5 

6 

.53 

.118 

.1414 

28 

k 

5 

6 

.53 

.148 

.14a 

29 

u 

5 

6 

.53 

.148 

.ail 

30-i;8 

C 

5 

6 

.1.8 

.148 
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METRIC  SYSTEM 


BASE  UNITS: 


Quantity 

Unit 

SI  Symbol 

Formula 

length 

metre 

m 

mass 

kilogram 

lig 

time 

second 

s 

electric  current 

ampere 

A 

thermodynamic  temperature 

kelvin 

K 

amount  of  substance 

mole 

mol 

luminous  intensity 

candela 

cd 

SUPPLEMENTARY  UNITS: 

plane  angle 

radian 

rad 

solid  angle 

steradian 

sr 

DERIVED  UNITS: 

Acceleration 

metre  per  second  squared 

mJs 

activity  (of  a radioactive  source) 

disintegration  per  second 

(disintegration)/s 

angular  acceleration 

radian  per  second  squared 

rad/s 

angular  velocity 

radian  per  second 

rad/s 

area 

square  metre 

m 

density 

kilogram  per  cubic  metre 

kg/m 

electric  capacitance 

farad 

K 

^■$rv 

electrical  condui:lance 

siemens 

S 

A/V 

electric  field  strength 

volt  per  metre 

V/m 

electric  inductanc  e 

henry 

II 

V-s/A 

electric  potential  difference 

volt 

V 

W/A 

electric,  reaiatanc  e 

ohm 

V/A 

electromoliva  force 

volt 

V 

W/A 

energy 

joule 

1 

N-m 

entropy 

loule  per  kelvin 

|/K 

force 

newton 

N 

kg-m/s 

frequency 

hertz 

Hz 

(cycle|/s 

illuminance 

lux 

lx 

Im/m 

luminance 

candela  per  square  metre 

cd/m 

luminous  flus 

lumen 

Im 

cd-sr 

magnetic  field  strength 

ampere  per  metre 

A/m 

magnetic:  flux 

weber 

Wb 

Vs 

magnetic  flux  density 

tesla 

T 

Wh/m 

magnetomotive  force 

ampere 

A 

power 

watt 

W 

|/* 

pressure 

pascal 

Pa 

N/m 

quantity  of  elec  tric  ity 

coulomb 

t: 

A-s 

quantity  of  heat 

loule 

1 

N-m 

radiant  interuity 

watt  per  steradian 

Wist 

specific  heat 

loule  per  kilogram-kelvin 

)/kg.K 

stress 

paacal 

Pa 

N/m 

thermal  conductivity 

watt  per  metre-kelvin 

Wim-k 

velcx.ity 

metre  par  aec;ond 

m/s 

visc.osity.  dynamic 

pascal-second 

Pa-s 

viscxMity.  kirtemalic 

square  metre  per  sec.ond 

m/s 

voltage 

volt 

V 

W/A 

volume 

cubic  metre 

m 

wavenumber 

reciprcxial  metre 

1 wava)/m 

work 

loule 

1 

N-m 

SI  PREFIXES: 


Multipliration  Kaclori 


I'rsfU  SI  Symbol 


1 000  (Kxj  rjoo  000 

lO'* 

I OTH)  000  000 

10* 

1 000  OOO 

10* 

I 000 

10* 

too 

10* 

10 

10' 

0 1 

10-' 

0.01 

10-* 

0 001 

10-* 

O OOO  001 

10-* 

0.000  000  001 

10-’ 

0.000  000  000  001 
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MISSION 
of 

Rome  Air  Develnjyment  Center 


RADC  plans  and  conducts  research/  exploratory  and  advanced 
development  programs  in  command/  control/  and  communications 
(C^)  activities/  and  in  the  areas  of  information  sciences 
and  intelligence.  The  principal  technical  mission  areas 
are  communications,  electromagnetic  guidance  and  control/ 
surveillance  of  ground  and  aerospace  objects/  intelligence 
data  collection  and  handling/  information  system  technology/ 
ionospheric  propagation/  solid  state  sciences,  microwave 
physics  and  electronic  reliability , maintainability  and 
compatibility . 
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